2006.13226v1 [astro-ph.HE] 23 Jun 2020 


e 
. 


arXiv 


MNRAS 000, 1-29 (2020) Preprint 25 June 2020 Compiled using MNRAS IATRX style file v3.0 


The Ly—Ly—Lradio relation and corona-disk-jet connection in 
optically selected radio-loud quasars 


S. F. Zhu,!?* W. N. Brandt,!?? B. Luo,^ Jianfeng Wu,’ 
Y. Q. Xue, 97 and G. Yang®? 


! Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 
? Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 
? Department of Physics, 104 Davey Lab, The Pennsylvania State University, University Park, PA 16802, USA 

4 School of Astronomy and Space Science, Nanjing University, Nanjing, Jiangsu 210046, China 

5 Department of Astronomy, Xiamen University, Xiamen, Fujian 361005, China 


6CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei £ 


7 School of Astronomy and Space Sciences, University of Science and Technology of China, Hefei 230026, China 
8 Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843-4242, USA 


? George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, Texas A&M University, College Station, TX 77843-42 


Accepted XXX. Received YYY; in original form ZZZ 


ABSTRACT 

Radio-loud quasars (RLQs) are more X-ray luminous than predicted by the X-ray- 
optical/UV relation (i.e. Lx cc L},) for radio-quiet quasars (RQQs). The excess X-ray 
emission depends on the radio-loudness parameter (R) and radio spectral slope (a;). 
We construct a uniform sample of 729 optically selected RLQs with high fractions 
of X-ray detections and a, measurements. We find that steep-spectrum radio quasars 
(SSRQs; o, < —0.5) follow a quantitatively similar Ly œ LY, relation as that for RQQs, 
suggesting a common coronal origin for the X-ray emission of both SSRQs and RQQs. 
However, the corresponding intercept of SSRQs is larger than that for RQQs and in- 
creases with R, suggesting a connection between the radio jets and the configuration 
of the accretion flow. Flat-spectrum radio quasars (FSRQs; a, > —0.5) are generally 
more X-ray luminous than SSRQs at given Luy and R, likely involving more physical 
processes. The emergent picture is different from that commonly assumed where the 
excess X-ray emission of RLQs is attributed to the jets. We thus perform model se- 
lection to compare critically these different interpretations, which prefers the coronal 
scenario with a corona-jet connection. A distinct jet component is likely important for 
only a small portion of FSRQs. The corona-jet, disk-corona, and disk-jet connections 
of RLGs are likely driven by independent physical processes. Furthermore, the corona- 
jet connection implies that small-scale processes in the vicinity of SMBHs, probably 
associated with the magnetic flux/topology instead of black-hole spin, are controlling 
the radio-loudness of quasars. 


Key words: quasars: general — X-rays: galaxies — galaxies: nuclei — galaxies: jets — 
black hole physics 


1 INTRODUCTION RLQs) are observationally distinguished from other quasars 
t d radio-quiet R b iri dio- 
Quasars are luminous active galactic nuclei (AGNs) whose jtaraed nadie aie ee USE Uy toute tod 

loudness parameter R = Lsanz/ L4499 à 2 10, where Ls5GHz 


central engines are supermassive black holes (SMBHs) that : A 
: à nan n and L are monochromatic luminosities at rest-frame 
are actively feeding. A significant minority of quasars (« 10— 


2096; e.g. Padovani et al. 2017) harbor a pair of power- 
ful relativistic jets that launch from a region close to the 
SMBH. Because relativistic jets are strong radio emitters, 
the quasars with powerful jets (termed radio-loud quasars, 


4400À 
5 GHz and 4400 A, respectively (Kellermann et al. 1989). 


The jets may be powered by the rotational energy of the 
SMBH and/or the inner accretion flow that is extracted by 
large-scale magnetic fields threading them (e.g. Blandford 
& Znajek 1977; Blandford & Payne 1982; see the review pa- 
per of Blandford et al. 2019). However, the mechanism that 
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triggers the production of powerful relativistic jets in only a 
minority of quasars is not clear. 

X-ray emission is nearly ubiquitous for quasars (Brandt 
& Alexander 2015, and references therein). The primary X- 
ray emission (~1-100 keV) from RQQs is radiated from a 
“coronal” structure containing hot plasma (e.g. Haardt & 
Maraschi 1993), which might be powered by magnetic recon- 
nection (e.g. Beloborodov 2017). UV photons from the inner 
accretion disk are up-scattered by electrons in the plasma 
to produce X-rays. This thermal Compton-scattering pro- 
cess leaves an imprint on the X-ray spectrum as a high- 
energy cutoff at —100-200 keV, which has been observed 
in local AGNs (e.g. Fabian et al. 2015) and a few high- 
redshift quasars (e.g. Lanzuisi et al. 2019). Interestingly, 
a non-linear correlation between the luminosities of the 
coronal structure and accretion disk has been established, 
Lov «* L5 A" which is referred to as the Loey—Ly so, 


relation (y = 0.6-0.8; e.g. Avni & Tananbaum 1986; Just 
et al. 2007; Lusso & Risaliti 2016).! The Ax—Ly soy A rela- 
tion describes the same non-linear correlation as a depen- 
dence of the shape of the optical/UV—X-ray spectral en- 
ergy distribution (SED) on disk luminosity. Here, aox = 
log(Lokev / L 5500 A) / log(vokev / 2500 A) is the two-point spectral 
index between rest-frame 2 keV and 2500 A (Tananbaum 
et al. 1979). 

The X-ray properties of RLQs are different from those 
of RQQs. RLQs are generally more X-ray luminous than 
RQQs of matched optical/UV luminosity (e.g. Zamorani 
et al. 1981; Worrall et al. 1987; Miller et al. 2011; Ballo 
et al. 2012). Their X-ray spectra are systematically flatter 
than those of RQQs (e.g. Wilkes & Elvis 1987; Reeves et al. 
1997; Page et al. 2005), particularly for the flat-spectrum ra- 
dio quasars (FSRQs; e.g. Grandi et al. 2006).? Furthermore, 
Compton-reflection features are weaker in RLQs (e.g. Reeves 
et al. 1997; Reeves & Turner 2000). Among low-redshift 
AGNSs, radio galaxies are more X-ray luminous than their 
radio-quiet counterparts (e.g. Gupta et al. 2018); broad-line 
radio galaxies (BLRGs) are found to have weaker reflection 
features than type-1 Seyfert galaxies (e.g. Wozniak et al. 
1998; Eracleous et al. 2000). However, no strong evidence 
supports the X-ray spectra of BLRGs being flatter than 
those of radio-quiet Seyfert galaxies (e.g. Sambruna et al. 
1999; Grandi et al. 2006; Gupta et al. 2018). The different 
X-ray properties of radio-loud and radio-quiet AGNs could 
be explained if the radio “core” of the jets contributes sig- 
nificantly in the X-rays a broadband component with a flat 
spectrum, in addition to the typical disk/corona emission be- 
ing present (e.g. Lawson & Turner 1997; Grandi & Palumbo 
2004).? However, except for a few FSRQs (e.g. Grandi et al. 
2006; Madsen et al. 2015), the X-ray spectra of BLRGs (e.g. 
Wozniak et al. 1998; Sambruna et al. 2009; Ronchini et al. 


! Throughout the paper, y with no subscript refers to the slope 
of the Ly—L,y relation for RQQs. 

2 The flat-spectrum and steep-spectrum objects have a, > -0.5 
and a, < —0.5, respectively. Here, a; is the power-law spectral 
index (i.e. f» cc v?*) in the radio band. 

3 The radio “core” here refers to the sub-arcsec component of the 
radio image that spatially coincides with the optical and X-ray 
(point-source) position of the quasar. It is presumably related to 
the base of the jet. 


2019) and steep-spectrum radio quasars (SSRQs; e.g. Lo- 
hfink et al. 2017) generally do not reveal a jet-linked flat 
continuum, suggesting that the orientation of the jets with 
respect to our line of sight might play an important role. 

The coupling between the coronal structure and the ac- 
cretion disk that is revealed by the Lokev—Ly soo À relation 
of RQQs probably exists for RLQs as well. Both the jet- 
launching region and corona are in the immediate vicinity 
of the SMBH, and connections (e.g. through a joint depen- 
dence on the magnetic field) between relativistic jets and 
the corona might be expected. The disks/coronae of quasars 
dissipate most of their radiated energy in the optical/UV 
and X-ray bands. Powerful relativistic jets have character- 
istic synchrotron radio emission, and might have contribu- 
tions in the X-rays. Therefore, X-ray, optical/UV, and radio 
are the key observational windows through which to peer at 
the nature of RLQs. An empirical Lokev —Lo5994 -L5Gnz rela- 
tion has long been sought for RLQs (e.g. Tananbaum et al. 
1983), in analogy to the Loxev—Ly sod relation for RQQs.* 
Indeed, the amount of the X-ray excess of RLQs over RQQs 
of comparable optical/UV luminosities increases with both 
radio-loudness parameter and radio luminosity (e.g. Miller 
et al. 2011), supporting the idea that the X-ray luminosity 
is determined by considering the power of both the disk and 
the jets, which are represented by L,<4)4 and Lsgg;, respec- 
tively. Except for extreme objects at extreme redshifts (i.e. 
quasars with log R > 2.5 at z 2 4; Wu et al. 2013; Zhu et al. 
2019), the Lokev-L» 599 4 ~L5GHz relation does not have an ap- 
parent redshift dependence (e.g. Worrall et al. 1987; Miller 
et al. 2011). 

However, studies of the Loxev 7L» 599 4 ~LSGHz relation for 
RLQs have generally lagged behind those of the Loyey—- 
Lso0À relation for RQQs (e.g. Worrall et al. 1987; Miller 
et al. 2011). On one hand, the sample size of RLQs is about 
an order of magnitude smaller than that of RQQs from opti- 
cal quasar surveys. On the other hand, to constrain the rela- 
tions to a comparable precision, studies of the Loxey—L 599 4 — 
Ls5gg; relation within at least a three-dimensional parame- 
ter space generally require the sample size to be larger than 
that of RQQs, for which a two-dimensional parameter space 
is sufficient. Furthermore, an extra dimension of RLQ prop- 
erties (i.e. radio luminosity) compared to RQQs also indi- 
cates an extra dimension of the model space. There are more 
candidate models for the Lokev ~La 599 4 -L5GHz relation from 
which to choose. Beaming effects of the jet emission might 
add another layer of complexity. 

The ambiguity in the functional form of the Loyey— 
L soo L5GHz relation makes its physical interpretation and 
implications unclear. It is possible that RLQs have a similar 
X-ray emitting disk/corona structure as that of RQQs and 
a distinct jet-linked X-ray component. The Loey—Ly 5994 — 
Ls5gg; relation would then simply describe the general pos- 
itive correlations of total X-ray luminosities with the opti- 
cal/UV and radio luminosities. 

Alternatively, the Lwev-Los994-Lscnuz relation might 
indicate a connection between the jets and disk/corona; e.g. 
the disks/coronae of quasars that have stronger relativistic 


^ Such empirical relations not only advance our understanding of 
AGN physics but also have broad practical applications in, e.g., 
SED fitting (e.g. Yang et al. 2020). 
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jets could be more X-ray luminous than those of quasars 
that have weaker or no relativistic jets. Such a connection 
might link AGNs to the phenomena of Galactic black-hole X- 
ray binaries (BHXRBs; or microquasars sometimes called) 
that are also powered by the black-hole accretion process 
and show couplings between jets and the accretion flow (e.g. 
Marscher et al. 2002; Merloni et al. 2003). Most BHXRBs 
are transients and show outbursts that last months to years 
(e.g. Remillard & McClintock 2006). During a typical out- 
burst, BHXRBs may cycle through (a few) state transitions 
that are marked by changes in spectral and timing prop- 
erties (e.g. Homan & Belloni 2005), as well as jet activity 
(e.g. Fender et al. 2004). The physical scales of SMBHs at 
the centers of massive galaxies are > 10° times those of the 
stellar-mass black holes of BHXRBs, making it practically 
difficult for observations to spot state transitions of indi- 
vidual AGNs directly (e.g. Schawinski et al. 2015). Instead, 
snapshot (as compared to the timescale of state transitions) 
observations across different wavelengths discover a great va- 
riety of AGNs (e.g. Padovani et al. 2017). While some of the 
varieties are caused by inclination-dependent geometry as we 
are able to observe only one aspect of each AGN (e.g. Netzer 
2015), accretion states of the central engine might also play 
an important role (e.g. Best & Heckman 2012). Investigating 
the disk/corona-jet connection of RLQs and establishing a 
phenomenological correspondence between AGN types and 
BHXRB states can shed light on the physics of black-hole 
accretion and relativistic jets. 

Previous works did not systematically compare these 
scenarios (e.g. Tananbaum et al. 1983; Zamorani 1983, 1984; 
Worrall et al. 1987; Miller et al. 2011). Specifically, they usu- 
ally focus on one functional model and obtain several sets of 
parameters that reflect different X-ray properties of different 
samples. Those sample-dependent empirical relations can be 
used to predict the X-ray luminosity for given optical/UV 
and radio luminosities, within a restricted parameter space. 
However, the driving mechanisms are hidden due to the lack 
of generality. Here we instead apply various models to the 
same sample and seek the most probable explanation of the 
data (i.e. model selection). Perhaps no single model is suit- 
able for all RLQs. For example, FSRQs and SSRQs might 
require separate mechanisms to explain their X-ray data. 
Then, we compare the results across all samples, investigat- 
ing their differences as well as similarities. 

We construct a large (> 700 objects) optically selected 
RLQ sample without regards to their radio/X-ray properties 
in § 2. Those RLQs span a broad parameter (i.e. luminosity, 
radio slope, and radio-loudness) space and have a high X-ray 
detection fraction (Z; 9096). Furthermore, almost all of them 
(Z; 96%) have basic radio spectral information. We perform 
model-independent, model-fitting, and model-selection anal- 
yses in 8 3. We compare with literature results and discuss 
physical implications in § 4. A summary of this paper and 
future prospects are in $ 5. In this paper, the quoted error 
bars represent lo uncertainties, and the upper limits are at 
a 9596 confidence level, unless otherwise stated. The spec- 
tral index (ay) follows the convention that fy œ v?^v. We use 
Lx, Luv, and Lradio interchangeably with Loyey, Lasoo» and 
Lsguz- The median statistic is widely used throughout the 
paper. We calculate medians using the Kaplan-Meier esti- 
mator (e.g. Kaplan & Meier 1958) in cases where the data 
contain non-detections. We use the bootstrapping method 
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if the uncertainties of medians are quoted. We adopt a flat- 
ACDM cosmology with Hy = 70 km s^! Mpc^! and Qn = 0.3. 


2 SAMPLE SELECTION 


We select new RLQs utilizing the Sloan Digital Sky Survey 
(SDSS; York et al. 2000). The radio data are from the Faint 
Images of the Radio Sky at Twenty-Centimeters (FIRST; 
Becker et al. 1995) and the NRAO VLA Sky Survey (NVSS; 
Condon et al. 1998). Archival Chandra and XMM-Newton 
observations are used to constrain the X-ray luminosities. 
The newly selected RLQs are combined with those from 
Miller et al. (2011) to form a final sample of 729 optically 
selected RLQs, which is summarized in Table 1. Compared 
with Miller et al. (2011), both the sample size and X-ray de- 
tection fraction are increased. Furthermore, we double the 
numbers of spectroscopic redshifts and radio slopes, the lat- 
ter of which affects the X-ray properties of RLQs (see 8 3). 
Importantly, the number of log R > 3 RLQs with reliable 
spectroscopic redshifts are significantly increased (by 70%). 
We will show in § 3.2.3 that such RLQs with the highest 
radio-loudness parameters have the largest statistical power 
in discriminating between models. 


2.1 New RLQs from the SDSS DR14 Quasar 
catalog 


2.1.1 Initial selection 


The SDSS DR14 Quasar catalog (DR14Q; Paris et al. 2018) 
covers a sky area of 9376 deg? and contains spectroscopically 
identified quasars from the Legacy Survey of SDSS-I/II, the 
Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS- 
III, and the extended Baryon Oscillation Spectroscopic Sur- 
vey (eBOSS) of SDSS-IV. The size of DR14Q (5.3 x 10° 
quasars) is a factor of ~ 7 times that of DR5Q (7.7 x 10+ 
quasars), which was utilized by Miller et al. (2011). The sky 
coverage of the FIRST survey has an extent of 10575 deg? 
and largely coincides with that of the SDSS. The NVSS cov- 
ers the entire sky north of Dec = —40 degrees but with a 
beam about 10 times larger than that of FIRST. Since the 
FIRST images have a better resolution, we select new RLQs 
from the matching results of DR14Q with the final catalog 
of the FIRST survey (Helfand et al. 2015). Considering the 
fact that the quasars in DR14Q are generally fainter than 
those in DR5Q, we only consider RLQs with log R > 2, which 
ensures mj < 21 quasars can be detected in the radio band 
given the (5c) flux limit of about 1 mJy of the FIRST sur- 
vey. The matching is performed as follows. We refer to each 
row in the FIRST catalog as a radio component. We adopt 
the method of Banfield et al. (2015) to distinguish resolved 
and unresolved components (cf. their Eq. 1 and Fig. 2). The 
radio flux of a resolved (unresolved) radio component refers 
to its integrated (peak) flux. We add the radio fluxes of all 
radio components within a radius of 90 arcsec around the 
optical position of a quasar and calculate a first correspond- 
ing radio-loudness parameter, which results in 24772 candi- 
date log R > 2 quasars in the redshift range of 0.5 « z < 4. 
Here, we have assumed o, = —0.5 to calculate Ls GHz from 
the observed 1.4 GHz flux. The i-band apparent magnitude 


(mi) is utilized to calculate L,,9, 4, where the K-correction 
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of Richards et al. (2006) and an optical spectral index of 
@ = —0.5 are assumed. 

In the above, a very large matching radius (90 arcsec) 
is adopted to ensure that the extended radio emission (e.g. 
from jets and lobes) associated with each quasar is recov- 
ered. However, in many cases, the R value calculated here 
is merely an upper limit, because background radio sources 
that are not associated with the quasar are also included. 
Visual inspection is required to eliminate such contamina- 
tion from background radio sources (e.g. Lu et al. 2007). 
To minimize the work of visual inspection, we first match 
the list of candidate log R 2 2 quasars with the observation 
catalogs of Chandra and XMM-Newton and apply unbiased 
empirical quality cuts. For the Chandra/ ACIS observations, 
we require 


Texp > 1000 + 35 x 10/8), (1) 


where Texp is the exposure time in seconds and @ is the off- 
axis angle of the quasar on the X-ray image, in units of 
arcmin. This criterion requires the exposure time to be at 
least 1 ks, and it requires additional exposure time at large 
off-axis angles to compensate for the loss of sensitivity due to 
the larger point spread function. Similarly, the quality cuts 
for XMM-Newton/EPIC-pn and XMM-Newton/EPIC-MOS 
are 


Texp > 1000 + 20 x 10/8)" (2) 
and 
Texp > 2000 + 20 x 10/8), (3) 


respectively. Note that 0 is bounded by the field of view 
of the telescopes (0 « 15’). In addition to the cut on the 
exposure time, we also require each quasar not to fall onto 
the edge of the detector or in CCD gaps. The requirements 
for sensitive X-ray coverage and quality cut here reduce the 
sample size to 1090. 

For these quasars, if there is no radio component in 
the annulus of 2 < r < 90 arcsec, no visual inspection is 
performed (330 quasars) and the association is assigned au- 
tomatically. We visually inspected the 4 x 4 arcmin? FIRST 
images of the remaining 760 quasars, where the radio com- 
ponents with apparent optical counterparts are labeled. Ra- 
dio components that are associated with other background 
sources in the field of view are discarded. In cases where real 
association exists (see Appendix A of Miller et al. 2011 for 
the utilized matching method), we associate the quasar with 
one or multiple radio components. 

Furthermore, even though FIRST has a deeper nom- 
inal flux limit than that of NVSS, it is known to be less 
sensitive to extended radio components (e.g. White et al. 
2007). Sometimes, the faint extended radio emission is “re- 
solved out" and completely missing in the FIRST catalog 
(e.g. Blundell 2003). The LOFAR Two-metre Sky Survey 
(LoTSS; Shimwell et al. 2017) aims to image the 120-168 
MHz northern sky with a sensitivity of 100uJy. The LoTSS 


5 By matching the DR14Q with FIRST, we might already miss 
RLQs with only faint and diffuse radio emission. However, those 
cases are rare (~ 2.3%; e.g. Lu et al. 2007) such that quasars with 
an extended morphology are almost always luminous in the radio 
band. 


data release 1 (DR1; Shimwell et al. 2019) covers a sky area 
of about 400 square degrees, which contains 63 of our RLQs. 
We compare the LoTSS and FIRST images of those RLQs 
and find three cases where FIRST resolves out extended ra- 
dio components, among which the FIRST fluxes are signifi- 
cantly lower than those of NVSS by 30%-50% for two RLQs. 
Even though those are rare cases (2/63 ~ 3%), we preferen- 
tially use the radio fluxes from NVSS over those of FIRST 
to avoid underestimating total radio emission, provided that 
the former is not contaminated by very nearby background 
radio sources. Note that we regard the FIRST component 
< 2 arcsec away from the quasar as the radio core. Using the 
ratio of the peak flux of the radio core to quasar total flux, 
we assess the dominance of the core (see Footnote 13). Given 
the angular resolution of the FIRST survey ( 5 arcsec), the 
core dominance here might overestimate the true core dom- 
inance that is revealed in very-long-baseline interferometric 
data. 

After the above procedures, the resulting sample has 
a size of 545. We label quasars as serendipitously observed 
if their X-ray observations have 0 > 1 arcmin, while the 
rest are labeled as targeted. We also label the quasars 
that are initially selected in SDSS color space, regard- 
less of their radio/X-ray properties; this allows us to con- 
struct an optically selected sample to match the method- 
ology of Miller et al. (2011). We consider those quasars 
spectroscopically confirmed in the Legacy Survey that have 
LEGACY TARGET! flags of “QSO”, “HIZ” and “serendip- 
itous” as optically selected, following Miller et al. (2011). 
For quasars targeted in BOSS, we consider the “CORE” and 
“BONUS” samples as indicated by the BOSS TARGETI 
flags. For quasars targeted in eBOSS, we include those 
belonging to the “CORE” sample as indicated by the 
EBOSS_TARGET1 flags. We refer readers to Paris et al. 
(2018) and references therein for the target selection flags of 
DR14Q. 


2.1.2 X-ray data analysis 


We first matched the quasar positions to the Chandra Source 
Catalog Release 2.0 (CSC 2.0; Evans et al. 2010) and the 
latest XMM-Newton Serendipitous Source Catalog (3XMM- 
DR8; Rosen et al. 2016) to obtain their X-ray fluxes (0.5-7 
keV for Chandra observations and 0.5-4.5 keV for XMM- 
Newton observations). Galactic-absorption correction is per- 
formed subsequently (Dickey & Lockman 1990). 

Then, for the observations that are not included in the 
source catalogs or where the quasar is not detected by the 
catalog pipeline, we download the observations and analyze 
the data manually. Specifically, data reduction, cleaning, and 
image extraction are performed using CIAO (v4.11) and SAS 
(v17.0.0) for Chandra and XMM-Newton observations, re- 
spectively. Raw source counts are extracted from a circular 
region centered at the quasar position and enclosing «9096 
of the total energy. Background counts are extracted from 
a source-free concentric annulus or nearby circular region. 


6 The target selection algorithms of the BOSS and eBOSS pro- 
grams might also utilize infrared photometry information where 
available. 
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Using the source and background counts of each observa- 
tion, we calculate the binomial no-source probability (Pg; 
Weisskopf et al. 2007), which is the conditional probabil- 
ity of producing counts that are equal to or larger than the 
observed counts in the source extraction region given the 
intensity of the background. We take a source as detected if 
Pg < 0.01 and calculate its net count rate. Instrumental re- 
sponse files (i.e. response matrix files and ancillary response 
files) at the position of the source are created using the stan- 
dard routines of CIAO and SAS; these include the appropriate 
aperture correction. With these instrumental response files, 
the net count rate is converted to X-ray flux using SHERPA. 
We have fixed the photon index to T = 1.5 for the power-law 
model, following Miller et al. (2011); the dependence of the 
X-ray flux on the choice of I is mild for reasonable choices of 
I. Note that the X-ray flux is corrected for Galactic absorp- 
tion by specifying the Ng value (Dickey & Lockman 1990) of 
the model in SHERPA. A source is treated as a non-detection 
if Pg > 0.01. We calculate the upper limits on the net counts 
at a 9596 confidence level using the algorithm of Kraft et al. 
(1991) for non-detected quasars, and the upper limits of the 
X-ray fluxes are calculated using the same method as detec- 
tions. 

The Galactic-absorption corrected X-ray flux (or upper 
limit), from either the source catalogs or manual analysis, is 
subsequently converted to the monochromatic luminosity at 
rest-frame 2 keV, Loyey. 


2.1.8 Radio spectral indices 


We matched the quasar list with sky surveys in the radio 
band to obtain radio slopes, using fluxes at 1.4 GHz and at 
another wavelength. Those radio surveys include the Green 
Bank 6-cm (GB6) Radio Source Catalog (4.85 GHz; Gre- 
gory et al. 1996), Westerbork Northern Sky Survey (325 
MHz; Rengelink et al. 1997), TGSS Alternative Data Release 
(150 MHz; Intema et al. 2017), and LoTSS DR1 (144 MHz; 
Shimwell et al. 2019). For quasars that are not matched with 
those radio surveys, we gather their multi-wavelength ra- 
dio fluxes in the NED’ and VizieR.? For the 78 quasars 
that do not have multi-wavelength radio data in the NED 
or VizieR, 60 of them are covered by the ongoing VLASS 
(Lacy et al. 2020), which provides the 3 GHz flux. We down- 
loaded VLASS quick-look images? of these 60 quasars and 
measured their 3 GHz fluxes using AEGEAN (Hancock et al. 
2012, 2018). Note that the peak and integrated fluxes from 
these quick-look images are corrected by factors of 1.15 and 
1.10, respectively (Lacy 2019). We inspect the FIRST/NVSS 
images while matching RLQs with the above radio surveys 
to ensure true association of radio components and to screen 
background contamination, which is important for the GB6 
data since it has a large beam size (> 3 arcmin).!° In total, 
263 quasars have high-frequency (» 1.4 GHz) radio data, 
while 346 quasars have low-frequency (« 1.4 GHz) radio 


7 https://ned.ipac.caltech.edu/ 

8 http://vizier.u-strasbg.fr/vizier/sed/ 

? https://archive-new.nrao.edu/vlass/quicklook/ 

10 We measure FIRST/NVSS-VLASS radio slopes for those 
quasars having FIRST/NVSS-GB6 measurements to ensure that 
the inhomogeneous beam sizes do not strongly bias the results. 
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data. Even though the radio slope measurements are not 
from homogeneous data sets, they are sufficient to distin- 
guish flat-spectrum and steep-spectrum objects. Further- 
more, for objects having fluxes at > 3 different frequencies, 
the radio spectral shape is generally consistent with a single 
power-law, which argues against a strong effect of redshift. 
However, we label those “GHz peaked” and related quasars 
(28 objects) that are inconsistent with a power-law spec- 
trum but seem to be peaked at GHz to 100 MHz frequencies 
(e.g. O’Dea 1998). Note that these objects are eventually 
excluded from the final sample since their X-ray properties 
are often different from typical RLQs (e.g. Siemiginowska 
et al. 2008). 


2.1.4 Cleaning the sample using multi-wavelength colors 
and optical/UV spectra 


We matched the 545 quasars resulting from initial selec- 
tion (8 2.1.1) with infrared sky surveys, including the Wide- 
field Infrared Survey Explorer (WISE; Wright et al. 2010), 
the UKIRT Infrared Deep Sky Survey (UKIDSS; Lawrence 
et al. 2007), the UKIRT Hemisphere Survey (UHS; Dye et al. 
2018), and the VISTA Hemisphere Survey (VHS; McMahon 
et al. 2013). We also matched our sample with the Pan- 
STARRS Stack Object Catalog (Chambers et al. 2016) to in- 
clude Y-band photometry. Including the SDSS, we have high- 
quality photometric coverage from the UV to mid-infrared. 
For the mid-infrared data, we have utilized the forced pho- 
tometry from the unWISE coadd images (Lang et al. 2016), 
which improves the number of detections in the W1-WA 
bands. Since not all the forced-photometry data are high- 
quality measurements, we replace those fluxes with signal- 
noise-ratio (SNR) smaller than 2 with their 9596 confidence 
upper limits, and flag them as non-detections.!! For those 
quasars that do not have forced photometry, we still use 
the measurements from the AIIWISE catalog. The forced 
photometry increases the number of detections by ~ 30—50 
for each mid-infrared band. These photometric data are cor- 
rected for Galactic extinction using the dust map of Schlafly 
& Finkbeiner (2011) before the following analysis. 

We found that some of the RLQs have abnormal multi- 
wavelength colors and SDSS spectra that are not consis- 
tent with those of typical quasars featuring big blue bump- 
dominated SEDs and strong optical/UV emission lines. We 
here further screen the sample based on their SDSS spec- 
tra and multi-wavelength colors. We matched the 545 RLQs 
with the SDSS-DR7 Quasar Catalog (DR7Q; Schneider et al. 
2010; Shen et al. 2011) and SDSS-DR12 Quasar Catalog 
(DR12Q; Paris et al. 2017). For the 501 quasars that are 
in DR7Q or DR12Q, 31 are classified as broad absorption 
line quasars (BALs). For the remaining 44 quasars, one is 
a BAL because its Balnicity index is larger than 0 (Páris 
et al. 2018). We have removed these BAL quasars from fur- 
ther consideration because their X-ray properties are likely 
affected by complex absorption (e.g. Miller et al. 2009). We 
also flagged 10 quasars whose SDSS spectra are Type II-like. 
They are relatively local (z < 0.76) with strong [O 11], weak 
Hf, and often red continua. Four of them are in the Type II 
quasar catalog of Reyes et al. (2008) or Yuan et al. (2016). 


11 http://wise2.ipac.caltech.edu/docs/release/allwise/expsup/sec3_la.html 
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Table 1. Summary of optically selected RLQs utilized in the paper. 


Sample No. of Sources X-ray Detections Serendipitous Spectroscopic z Radio Slope 
All RLQs 729 657 (90.1%) 622 (85.3%) 587 (80.5%) 704 (96.6%) 
FSRQs 394 363 (92.1%) 340 (86.3%) 319 (81.0%) 394 (100%) 
SSRQs 310 275 (88.7%) 258 (83.2%) 253 (81.6%) 310 (100%) 


Table 2. Properties of optically selected RLQs utilized in the paper, in ascending order of RA. Only the top 5 objects are listed. 


SDSS Name z mi logLyso)4 logLscuz logLkey logR XFlag* sdP spec? Qr cD4 
000442.18+000023.3 1.008 18.98 30.22 32.09 26.91 1.74 1 1 1 - 0.59 
000622.60—000424.4 1.038 19.58 30.05 34.94 27.32 4.77 1 1 1 -0.55 0.80 
001646.54—005151.7 2.243 21.00 30.20 32.66 26.36 2.34 1 1 1 -0.21 1.00 
001910.95+034844.6 2.022 20.30 30.35 32.91 26.54 2.44 1 1 1 —0.43 0.56 
003054.63+045908.4 2.201 20.95 30.19 33.81 26.59 3.50 1 1 1 -0.41 0.76 

? Tf the quasar is detected in X-rays, XFlag = 1, while XFlag = 0 if otherwise. 

P sd = 1 if the X-ray observation is labeled serendipitous and sd = 0 if the quasar was targeted. 

* If the quasar is spectroscopically confirmed, spec = 1. If the redshift is based on photometric data, spec = 0. 

4 The ratio of the flux of the radio core over the total radio flux (see § 2.1.1). 

These quasars also likely have complex X-ray absorption, so Ass (À) 
105 104 103 


are screened. We screened out quasars that have red SDSS 
and/or infrared colors that satisfy the following criteria: 


u — W4 > 12.0 or u— W3 > 9.4 or u — i > 0.95, z < 1.5 
(4) 

g- W4 > 11.9 or g—W3>9.7 or g — i > 1.0, 1.5<z<2.5 
(5) 
r-W3»84orr-z»0.6 2.5«z x4 (6) 


Note that these color cuts effectively remove red quasars (79 
objects) that either suffer from severe dust extinction or have 
prominent jet emission in the infrared through UV bands. 
Furthermore, we exclude quasars that do not have strong 
emission lines, since their rest-frame optical/UV emission 
might be contaminated by strongly boosted jet emission, 
rendering the observed L,:,, 4 an overestimate of their disk 
power. Specifically, we measure the rest-frame equivalent 
width (REW) of Mg 1 (when it is covered by SDSS spectra) 
and exclude those quasars (3 objects) with REW < 10 A. Af- 
ter these procedures, the resulting sample has a size of 421, 
among which 327 quasars are selected using SDSS colors 
(see the last paragraph of § 2.1.1). Note that serendipitous 
X-ray coverage makes up about 79% of this sample, so the 
majority of these quasars have X-ray observations unbiased 
by source properties. Those quasar jets with the most ex- 
tremely boosted emission are under-represented in our sam- 
ple, due to the color cuts and the constraint on emission-line 
strength. However, their fraction is estimated to be only 
S 6%, considering a conservative Lorentz factor of l = 8 
for the jets and a half-opening angle of 30 deg for the dusty 
torus, and thus the results of this paper will not be strongly 
biased. 


2.2 Improving the Miller et al. (2011) RLQ 
sample 


Miller et al. (2011) construct a sample of 791 quasars with 
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Figure 1. The composite median SED (orange) of optically se- 
lected RLQs, from the mid-infrared to UV. The top and bot- 
tom ticks are in units of rest-frame frequency and wavelength, re- 
spectively. The small dots are data points for individual quasars, 
where grey symbols are non-detections in the WISE bands, distin- 
guished from detections (black). Any data points with rest-frame 
frequency > vj5,44 are not used in constructing the composite 
SED due to intergalactic absorption. Also plotted is the compos- 
ite SED of the local RLQs from Shang et al. (2011, blue) for 
comparison. 


log R* > 1.!? The overall X-ray detection fraction of this 
sample is high (8596). Their full sample consists of 654 op- 
tically selected primary RLGs, and the rest are supplemen- 
tary RLQs that are generally not optically selected. The 


12 Miller et al. (2011) use R* = Lscuz/Lys9)4 to quantify radio- 
loudness, which can be converted to R using log R = log R* —0.123. 
In Miller et al. (2011), quasars with 1 < log R* < 2 and log R* > 2 
are referred to as radio-intermediate quasars (RIQs) and RLQs, 
respectively. We call all their quasars RLQs, even though some of 
them have log R values slightly smaller than unity. 
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Figure 2. The composite median SDSS spectra (in Afa represen- 
tation) of optically selected RLQs. From top to bottom, RLQs 
are divided into bins of redshift, radio slope, and radio-loudness. 
The dash-dotted vertical lines show the frequencies of six emis- 
sion lines as labeled. Note that the composite spectra are arbi- 
trarily shifted vertically to avoid overlapping. In all cases, promi- 
nent emission lines are apparent, indicating that the optical/UV 
emission of the RLQs is not contaminated by strong boosted jet 
emission. 


X-ray data for the primary sample are from archival Chan- 
dra, XMM-Newton, and ROSAT observations, most (86%) 
of which are serendipitous and thus unbiased with respect 
to RLQ properties. 

However, about half of their primary RLQs lack radio 
spectral information. Furthermore, the primary sample of 
Miller et al. (2011) includes 312 spectroscopically confirmed 
quasars culled from the SDSS DR5 Quasar catalog (DR5Q; 
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Schneider et al. 2007), while the remaining 342 quasars are 
photometrically selected but not spectroscopically confirmed 
(Richards et al. 2009). The improvement upon RLQs from 
Miller et al. (2011) thus mainly comes from higher fractions 
of radio-slope measurements and spectroscopic redshifts. 

We measure the radio slopes using the same method as 
for newly selected RLQs by matching with radio sky surveys, 
NED and VizieR, and the quick-look images of VLASS. The 
spectral indexes are measured for 719 quasars, making up 
91.1% of all quasars from Miller et al. (2011). As to the 
optically selected RLQs, their radio-slope measurements are 
91.4% complete. Note that to obtain a final data set with 
uniform measurements in the radio bands, we also update 
the total radio fluxes for RLQs from Miller et al. (2011) 
using NVSS as in § 2.1.1. 

For the 342 photometrically selected quasars in the pri- 
mary sample of Miller et al. (2011), SDSS/BOSS/eBOSS 
spectra from the SDSS/BOSS spectrograph can now be 
found for 183 of them. We exclude 9 quasars whose rest- 
frame optical/UV spectra are BL Lac-like with a feature- 
less continuum. Another 9 quasars are excluded because 
they are likely to be BAL quasars. For the rest (165) of 
the quasars with new spectroscopic measurements, we re- 
place their photometric redshifts with spectroscopic red- 
shifts and update their luminosities accordingly. Note that 
catastrophic failures in the photometric redshifts where 
[Zphoto — Zspec|/(1 +Zspec) > 0.15 make up x5% of these quasars, 
suggesting that the photometric redshifts of the remaining 
(159) quasars that still lack spectroscopic measurements are 
largely reliable. 

As for the new RLQs, we also gather multi-wavelength 
data (mid-infrared to UV) for the sample of Miller et al. 
(2011) and apply color cuts described in § 2.1.4, which elim- 
inate 93 objects. Another 9 objects are excluded since their 
X-ray observations do not satisfy the criteria in § 2.1.1. 


2.3 The final optically selected RLQ sample 


We merge the newly selected sample with the refined sam- 
ple of Miller et al. (2011), resulting in 940 unique RLQs. 
Note that some quasars might be present in both samples, 
mainly because many photometric quasars of Miller et al. 
(2011) now have optical/UV spectra. For the physical quan- 
tities describing these quasars, we use their newly measured 
values in this paper. Since some quasars in the SDSS QSO 
catalog are selected by their radio or X-ray properties, we 
utilize only the optically selected subset (789 objects) that 
does not suffer from apparent radio/X-ray selection biases 
in the following sections. Note that 49 quasars are not con- 
sidered since they either have peaked radio spectra or are 
compact objects with steep spectra. P? Another 11 quasars 
are excluded as rare cases that are known to be X-ray weak 
(2 quasars; e.g. Risaliti 2005), are in the vicinity of luminous 
clusters (1 quasar; e.g. Rykoff 2009), or have exceedingly flat 
X-ray spectra signifying strong absorption (8 quasars; e.g. 


13 Compact steep spectrum objects generally show a radio spec- 


trum that is peaked around 100 MHz. We exclude those quasars 
with core dominance > 0.9 and radio slope < —0.8, even if a peak 
is not revealed. 
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Figure 3. The optically selected RLQ sample in the Lyagio-Luy plane, where the data points are color-coded by their redshifts. Different 
symbols have been used to indicate the radio spectral information. Five lines defined by constant radio-loudness parameters are also shown. 
We also show 59 RQQs from Laor & Behar (2008), where squares and downward arrows are detections (50 quasars) and non-detections 


(9 quasars) in the radio band, respectively. 


Miyaji et al. 2006). The properties of the resulting 729 RLQs 
are listed in Table 2. 


We show in Fig. 1 the infrared-to-UV composite SED of 
the final optically selected RLQs, which has the blue color 
of typical Type I quasars that are dominated by the emis- 
sion from the accretion disk. Note that the composite SED 
(orange) in Fig. 1 utilizes ~ 10^ photometric data points of 
» 700 quasars. In comparison, we show the composite SED of 
RLQs from Shang et al. (2011) in the same plot. Fig. 2 shows 
the composite median spectra for various subsets of optically 
selected RLQs. They are divided into redshift, radio slope, 
and R bins in the three panels of Fig. 2. Prominent emission 
lines are present in all composite median spectra, indicating 
that the optical/UV continuum of the RLQs is generally not 
contaminated by boosted jet emission. Indeed, from the cat- 
alogs of Shen et al. (2011) and Paris et al. (2017), 573/729 
quasars of the final sample have REW measurements for at 
least one emission line among (C Iv, Mg ri, H£), and the me- 
dian REWs are (4253, 372, 77*6) A and (48*3,40*1,83*8) A 
for FSRQs and SSRGs, respectively. This indicates that even 
the FSRQs in our sample do not have a substantial opti- 
cal/UV continuum contribution from a jet. 


We show the final sample in the Lsggz-L plane in 


2500À 


Fig. 3. The color of each RLQ indicates its redshift, accord- 
ing to the color bar on the right-hand side. Furthermore, we 
show four lines on the plane that are defined by constant 
radio-loudness parameters. 


2.4 Comparison optically selected RQQ samples 


To compare the X-ray properties of RLQs with those of 
RQQs, we also utilize a large sample of RQQs from Lusso 
& Risaliti (2016). These RQQs were selected from the 
SDSS quasar catalog, and the X-ray data are exclusively 
from XMM-Newton observations. Following Lusso & Risal- 
iti (2016), we select RQQs from their Table 1 and Table 2 
that satisfy E(B-V) < 0.1, S/N > 5, and 1.9 < Tx < 2.8 from 
the main sample (cf. Table 3 of Lusso & Risaliti 2016). We 
further exclude a small number of RQQs that are not opti- 
cally selected. The resulting sample has a size of 1074, among 
which 699 have detected X-ray emission. We use this sam- 
ple to constrain the Ly—Lyy relation, which is established for 
RQQs that span > 4 decades in luminosity (e.g. Just et al. 
2007; Lusso & Risaliti 2016) and show no redshift depen- 
dence up to z > 6 (e.g. Vito et al. 2019). 

Furthermore, we utilize z « 0.5 PG quasars from 
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Figure 4. The RQQ samples in the Ly—L,, plane. The black 
symbols are (SDSS) RQQs from Lusso & Risaliti (2016), while 
the (PG) RQQs from Laor & Behar (2008) are further divided 
into two bins (orange squares and red diamonds) by their median 
radio-loudness parameter (R = 0.25). Fitting a line to the black 
data points results in the equation at the lower-right corner (see 
Table 4). The uncertainties of this relation are indicated by the 
black-shaded region. The PG quasars are consistent with the fit- 
ted Lx—Luy relation. No apparent radio dependence for PG RQQs 
is found. 


Boroson & Green (1992), which are optically selected and 
have relatively deep radio constraints (e.g. Kellermann et al. 
1989, 1994). We only consider a subsample of 59 RQQs 
that do not have strong Civ absorption from Laor & Behar 
(2008), 50 of which have detected radio emission while the 
remaining 9 have upper limits. All 59 quasars are detected 
by ROSAT in X-rays (Brandt et al. 2000; Steffen et al. 2006) 
except for one (i.e. PG 1259+593), which is instead detected 
by XMM-Newton (Ballo et al. 2012). 

The two RQQ samples are shown in Fig. 4, where the 
RQQs from Lusso & Risaliti (2016) are shown as black dots 
(detections) and downward arrows (non-detections). We fit 
the Lx—Luy relation of this sample using the method de- 
scribed in § 3.2.2. The results are given in both Fig. 4 and 
Table 4, which are consistent with those of Lusso & Risaliti 
(2016). The 59 PG quasars are consistent with the Lx—Luy 
relation in Fig. 4. The 20th and 80th percentiles of their 
radio-loudness parameters are R = 0.10 and R = 0.64, respec- 
tively. We divide these PG quasars into two bins separated 
by their median radio-loudness parameter (R = 0.25) to as- 
sess potential radio dependence, which is not found in Fig. 4. 
However, we do not treat this statement as conclusive since 
the sample has a small size and contains only low-luminosity 
quasars. 


3 THE RELATION BETWEEN X-RAY, 
OPTICAL/UV, AND RADIO LUMINOSITIES 


3.1 Insights from scatter plots 


In the top panel of Fig. 5, the Ljijgo-Luy plane is divided 
into 3 x3 sub-regions, the boundaries between which are the 
33rd and 66th percentiles along each axis. For each bin, the 
median Lx is calculated and indicated by the depth of the 
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Figure 5. Top: The dependence of X-ray luminosity on log Lsauz 
(y-axis) and log Lo, (x-axis). The scattered points are all op- 
tically selected RLQs, which are further grouped into nine bins. 
Those bins are separated by the 33rd and 66th percentiles along 
each axis. We calculate the median X-ray luminosity of each bin 
as indicated by the color bar on the right-hand side. Clearly, Ly 
increases with both increasing Lsgy, and L soo- Bottom: Same 
as the top panel for the R-L,,,,; plane. 


È FSRQs | 
B^ SSRQs 


Figure 6. The X-ray luminosities of FSRQs (green dots) and 
SSRQs (purple squares) over those of RQQs at given Ly soo are 
plotted against log R. We group each sample into three R bins 
and calculate their medians and uncertainties. 


color according to the color bar on the right-hand side. The 
median Lx increases in the positive direction for both Lyagio 
and Luy. We repeat the analysis for the R-Luv plane in Fig. 5 
(bottom), where the distribution of data points is more ho- 
mogeneous than that in the upper panel of Fig. 5. Both 
figures indicate that there are strong radio and optical/UV 
dependences for the X-ray luminosities of RLQs. 

We investigate the radio dependence for FSRQs and 
SSRQ separately in Fig. 6, where their X-ray luminosities 
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Figure 7. The L,—L,, relations for all RLQs, FSRQs, and SSRQs, 
from top to bottom. The resulting slope with uncertainty estima- 
tion is at the lower-right corner of each panel. The Lx—Lyy relation 
for RQQs (black) is also plotted for comparison. At given opti- 
cal/UV luminosity, RLQs are more X-ray luminous than RQQs 
in all panels. The slope of SSRQs is consistent with that of RQQs 
(y = 0.63*0-07). while the slope of FSRQs is steeper at a ~ 2.50 


y= -0.02 
significance level. 


divided by those of RQQs at given Ly, (using the Ly—Lyy 
relation in Fig. 4) are shown as a function of log R. FSRQs 
and SSRQs are each divided into three R bins by the 33rd 
and 66th percentiles of all RLQs (log R = 2.25 and 2.78). The 
medians with estimated uncertainties for each bin are calcu- 
lated. In Fig. 6, the X-ray luminosities of both RLQ types 
are larger than those of RQQs and show strong dependence 
on log R. The most X-ray luminous objects at each log R 
are almost always FSRQs, probably due to their seemingly 
larger scatter than that of SSRQs. Furthermore, the medians 
of FSRQs are larger than those of SSRQs by 1096-2096. 
The Lx—Luy relations for RLQs and RQQs are compared 
in Fig. 7 (top). The slope of the regression line for RLQs is 


slightly steeper than that for RQQs (i.e. y = 0.63+0-02), Note 


that the fitting method here is the same as for Fig. 4 and 
is described in § 3.2.2. In the middle and bottom panels 
of Fig. 7, FSRQs and SSRQs are compared with RQQs, re- 
spectively. The Lx—Lyy relation for FSRQs (slope= 0.74+9-04) 
appears steeper (at about 2.50 significance) than that for 
RQQs, while the slope for SSRQs (slope= 0.61+9-04) is con- 
sistent with y (to within ~ 3%). Therefore, the X-ray lumi- 
nosities of FSRQs and SSRQs probably have different radio 
and optical/UV dependences (at about 2.30 significance) as 
indicated by Fig. 6 and Fig. 7. 

In short, the X-ray luminosity of RLQs depends on Lyy, 
Lyadio (or R), and radio slope,!* which has also been con- 
cluded by many previous works (e.g. Worrall et al. 1987; 
Brinkmann et al. 1997; Grandi et al. 2006; Miller et al. 2011). 
In addition to the disk-corona interplay that is revealed by 
the Ly—-L,y relation for RQQs, we probably need at least 
two more mechanisms, or one mechanism that is affected by 
two key parameters, for those dependences in RLQs. We fo- 
cus on the Lx—Lyy relation and consider R and radio slope 
as controlled additional factors. Specifically, we fit the Lx- 
Luy relations for the three R bins of FSRQs and SSRQs (as 
in Fig. 6), separately, and show the results in Fig. 8. The 
corresponding 1c confidence regions are in Fig. 9. 

A notable pattern for SSRQs emerges that the slope of 
the Lx—Luy relation is always consistent with that for RQQs, 
while the intercept increases monotonically with R. If the ex- 
cess X-ray luminosity of SSRQs relative to RQQs is caused 
by emission from the core region of the jets, the contribution 
of this component increases with R, as expected. However, 
two properties are required further by Fig. 8 (bottom) that 
(a) the jet component depends strongly on Luy and (b) the 
slope of this dependence is consistent with y. For example, 
the SSRQs in the most radio-loud bin are on average a fac- 
tor of about 3 more X-ray luminous than RQQs (see Fig. 6), 
which means that the jet component (if it exists) dominates 
their X-ray luminosities; the dependence of their Lx on Luy 
is still strong with a slope = 0.60+0-05 (see Fig. 9), which is 
only different from y by ~ 5%. Indeed, perhaps we do not 
need two distinct X-ray components (i.e. from the corona 
and jet core) that are indistinguishable with regard to their 
correlations with Luy and their spatial properties.!5 A more 
natural and compact explanation is that the X-ray emission 
of SSRQs reveals only one component and is produced in 
the same way as for RQQs.; this component is attributed to 
the hot corona that is coupled with the disk (i.e. the Lx— 
Luay relation). Furthermore, there is a connection between 
the activity of the corona (i.e. the intercept of the Lxy—Lyy 
relation) in the immediate vicinity of the SMBH and that 
of the large-scale jets (i.e. R), such that those RLQs harbor- 
ing more powerful relativistic jets also have more powerful 


coronae. 16 


1^ The radio slope reflects, roughly, the direction of the jet relative 
to our line of sight. 

15 The current X-ray telescope with superb sub-arcsec resolution 
(i.e. Chandra) cannot resolve the nuclear X-ray emission of RLQs 
into different spatial components. 

16 In principle, another possible explanation for Fig. 8 (bottom) 
could be that even RQQs have a jet-linked X-ray continuum (as 
many RQQs do have small jets, and some envision the X-ray 
corona in radio-quiet objects to be the base of a jet). Following 
this explanation, we should expect a corona-jet connection for 
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Figure 8. The L,—L,, relation for FSRQs (top) and SSRQs (bottom) divided into three R bins as indicated by the color-bar on the 
right-hand side. The slopes and intercepts for the relations are given in Fig. 9. The Lx—Luy relation for RQQs (black) is plotted for 
comparison. For SSRQs, the slope is always consistent with that of RQQs while the intercept increases with R. FSRQs follow another 


pattern where both the intercept and slope increase with R. 


'The results for FSRQs are not as clear in Fig. 8 and 
Fig. 9. For the lowest R bin, the slope is consistent with 
y, supporting the idea that the observed X-rays are also 
dominated by the corona, and the intercept is consistent 
with that of SSRQs within error bars. However, the slope 
increases with R to 0.7-0.9 for FSRQs in the other two 
bins, and the intercepts are also larger than for the corre- 
sponding SSRQs (although only suggestively for the second 
R bin). Therefore, the Lx—Luy relation for the moderate-to- 
most radio-loud FSRQs probably requires additional mech- 
anisms. However, we stress that the medians of FSRQs and 
SSRQs of similar radio-loudness in Fig. 6 are only different 
at a 10%-20% level. Therefore, the additional mechanisms 
affecting FSRQs are likely secondary such that their X-ray 


RQQs as well, which is, however, not supported (see § 3.3 and 
Fig. 14). More generally, luminosity correlations and X-ray spec- 
tral properties support a thermal Compton-scattering origin for 
the X-ray continuum of RQQs (see § 1). 
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luminosities are generally controlled by the enhancement of 
the putative coronal emission as for SSRQs.!7 We test in the 
next section whether a distinct component from the jet core 
might play an important role or not. 


3.2 Parameterized modeling 


Qualitatively similar plots to those in § 3.1 can also be found 
in previous works (e.g. Worrall et al. 1987; Miller et al. 2011), 
although our quantitative constraints here are considerably 
tighter owing to our improved data quality. For example, the 
results in Fig. 9 are in line with those in Fig. 1 of Worrall 
et al. (1987). However, despite such consistent plots, this 
paper and the literature reach different conclusions. Specifi- 
cally, we suggest that the disk/corona of RLQs are systemat- 
ically more X-ray luminous than those of RQQs, while both 


17 In consequence, this statement is probably true for almost all 


RLQs. 
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'Table 3. Models used to explain the X-ray luminosity of RLQs. 


Functional Model 


m Tradi 
Decomposition? : L, = ALIY + 9 [radio 


radio 

A Tuv 8 l'radio 
Radio-loud quasars 
I: log Lo kev = log Lysoq à + adio 108 L 102 (2509 Jade? pay 4 - 

: 108 E2 kev = Œ + Yuv 108 L2500 A + Yradio 08 45 GHz 4400 radio Yuv + Yndio 
II: log Ly pey = log (AL gra ) A N Yav B Yradio 
Yra io uv radi. Yradio ^" A 2 

II: log Lo yv = log (AL; odo 12. + BL? atte) a(g) "^ "mae yav thao B Yadio 
Radio-quiet quasars 
log Laxey = B + y log Ly soa 108 y - - 


? The X-ray emission from RLQs is divided into a component that is from the disk/corona (ALI) and 


another component that is from the jet core (BL, 


Tradio 
adio 


), the latter of which may not be explicitly present 


in all the models (e.g. Model I). The quantities that are supposed to have the same physical meaning 
and thus be comparable between different functional models are also listed. 
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Figure 9. The lo confidence regions for slope and intercept of 
the Ly—Lyy relations in Fig. 8. 


Worrall et al. (1987) and Miller et al. (2011) attribute the 
excess X-ray emission of RLQs to the jet core. These two in- 
terpretations need to be further compared. In this section we 
use formal model selection to show that our interpretation 
in § 3.1 is indeed preferred by the data. Furthermore, the 
varying contributions of the two additive X-ray components 
(jets and corona) are probably more easily revealed by direct 
model fitting than by finding the proportionality of total X- 
ray luminosity with luminosities at other wavelengths (e.g. 
§ 3.1). 


3.2.1 Models for the X-ray-optical/U V-radio relation of 


RLQs 


We first set up a decomposition for the X-ray emission from 
RLQs, and then discuss three functional models within this 
framework. We write the X-ray luminosity of RLQs as 
Lx = ALY + Br nds 


radio * 


(7) 
where ALY and Brom represent the parts of the emis- 
sion from the corona and jets, respectively. Here, the X-ray 
luminosity from the disk/corona depends on the disk lumi- 


nosity (Luv) to the power of [yy and is subject to a normal- 


ization factor A. This Ly-L,y relation for the disk/corona 
emission of RLQs has thus been assumed to be of similar 
functional form to that of RQQs. Whether the specific pa- 
rameters of this relation are universal to both RLQs and 
RQQs is not assumed. Given that the X-ray and radio lumi- 
nosities of the sample span 4-6 decades, we adopt a power- 
law dependence for the case of the jet component as well (i.e. 
Lx jet = Ide This functional form is reasonably flexible 
and can accommodate many potential underlying physical 
emission mechanisms. In practice, the data might not re- 
quire both components to be present, in cases where one of 
the normalization factors is consistent with zero (most likely 
B). 

In Table 3, we list three functional models that describe 
different Lev -Lasog A Lsauz relations, all of which are in- 
terpreted in the context of Eq. 7. For example, when we 
assess the dependence of the X-ray luminosity on the opti- 
cal/UV luminosity of the disk/corona, i.e. Tuy, we will use 
Yuv + Ys dio of Model I and Model III, while we will use yuy 
of Model II. 

Model I was utilized by previous works because it de- 
scribes a joint dependence of the X-ray luminosity on both 
the optical/UV and radio luminosities (see Fig. 5) with 
the smallest number of parameters (e.g. Tananbaum et al. 
1983; Worrall et al. 1987; Miller et al. 2011). To preserve 
consistency with past work, we still use the form Loyey « 

Yuv . Lodo 
2500À SGHz ; 
using another form that Loyey « RY% aio Lado 18 We note 


from Fig. 8 that, for fixed R, the i.i. relation of RLQs 
has a slope such that Tuy = yuy + Via dio ^ Y for the majority 
of bins. Previous model fitting results using Model I (e.g. 
Table 7 of Miller et al. 2011) are roughly consistent with 
this suggestion. Here, a corona-jet connection is parameter- 
ized so that the normalization factor (intercept) correlates 


with the radio-loudness parameter to the power of Via dio (i.e. 


in model fitting, while we interpret the results 


18 Another advantage of this preference in model fitting is that 
there is non-zero covariance between the measurement errors of 
R and L soo» while Lsgg; and L are independent measure- 
ments. 


2500. 2500À 
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A x RY ratio). An X-ray emitting region that is associated with 
the relativistic jets is not explicitly included (i.e. B is set to 
zero). 

Similar to Eq. 7, Model II explicitly divides the X-ray 


emission of RLQs into corona (AL oA) and jet core 


(BI, components. However, it stands for a special but 


commonly assumed scenario where ADS x might be con- 
sistent with the X-ray luminosity of RQQs; therefore, RLQs 
are treated as a basic combination of a radio-quiet quasar 
engine with additional powerful jets from the perspective of 
the X-rays (e.g. Worrall et al. 1987). For the jet component, 
the X-rays might correlate linearly (i.e. Yradio = 1) with the 
radio emission from the same region as suggested by previous 
works (e.g. Zamorani 1984; Browne & Murphy 1987; Wor- 
rall et al. 1987). No apparent connection between the jets 
and the corona is present in this model; each contribution 
is mathematically independent. Model II has been applied 
to FSRQs (e.g. Zamorani 1984; Worrall et al. 1987) and SS- 
RQs (e.g. Zamorani 1984) in previous works. With small-size 
samples, their results actually hinted at a difference between 
the coronae of RLQs and those of RQQs, though no solid 
conclusion was reached.!? 

It is also possible that both Model I and Model II are 
describing part of reality. In 8 3.1, we require a corona-jet 
connection for the differences between RQQs and SSRQs and 
cannot rule out a jet component in FSRQs. We thus propose 
Model III (see Table 3) that combines those features. Model I 
and Model II are special cases of Model III, where either B 
Or o sm is set to zero, respectively. From this point of view, 
we utilize only Model III and discuss whether its certain 
parameters are zero or not below. 


3.2.2 Methods 


Following Miller et al. (2011), we first normalize the quasar 
luminosities to log Lasoo 4 = 30.5, log Ls Guz = 33.3, and 
log Lo key = 27.0, which are near to the median luminosities 
of our sample. The maximum-likelihood estimates are then 
obtained for the models in Table 3. The model fitting is per- 
formed using LMFIT (Newville et al. 2014), and the likelihood 
function takes into account the upper limits on the X-ray lu- 
minosities (e.g. Isobe et al. 1986). The maximum-likelihood 
and uncertainty estimations of model parameters are cal- 
culated using a Markov chain Monte Carlo (MCMC) algo- 
rithm (EMCEE; Foreman-Mackey et al. 2013). This method 
is thus equivalent to a Bayesian approach with flat priors. 
We initially leave all parameters free to vary to reveal the 
most general results. However, we fix the parameter values 
for cases where the parameters are either unconstrained or 
unphysical (see details in Appendix A), which also prevents 
the model performance from being imprecisely characterized 
by the model-selection methods we use below (e.g. Nelson 
et al. 2020). 

Second, as to model selection, the likelihood-ratio test 
(e.g. F-test) is not applicable here since Model I and Model II 
are not nested to each other (e.g. Protassov et al. 2002). We 


1? Zamorani (1984) does not provide error bars. Worrall et al. 
(1987) have a slightly larger sample than that of Zamorani (1984), 
but the constraint is still statistically insignificant. 
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Figure 10. The best-fit results of all optically selected RLQs 
using Model I (blue) and Model II (orange). Model I defines a 
plane while Model II defines a curved surface. 


thus use standard information criteria (ICs) to compare the 
performance of models. We adopt the widely used Akaike 
Information Criterion (AIC; Akaike 1974) and Bayesian In- 
formation Criterion (BIC; Schwarz 1978), which are defined 
as 


AIC = -21n Lmax + 2k, (8) 
BIC = -21n Lmax + kIn N. (9) 


Here, Lmax is the maximum likelihood, k the number of free 
parameters of the model, and N the number of data points. 
Under such definitions, the model with the smallest infor- 
mation criterion is selected. ICs thus favor high model like- 
lihood (Lmax) and penalize high model complexity (k). BIC 
penalizes model complexity more severely than AIC when 
N > 7. The significance of the selection result (i.e. the ev- 
idence for the selected model) is indicated by AAIC (e.g. 
§ 2.6 of Burnham & Anderson 2002) and ABIC (e.g. § 3.2 of 
Kass & Raftery 1995) between models. In this paper, AIC 
and BIC always agree on the best model, and in most cases, 
the best models are strongly favored with |AAIC| = 5 and 
|ABIC| = 5. It is probably not safe to rely solely on statistical 
methods and to draw conclusions without considering phys- 
ical plausibility. The results of the following two subsections 
are secured by the fact that the ICs-selected model fits also 
depict the most reasonable physical pictures. 


3.2.9 Results of all RLQs without distinguishing radio 
slope 


The fitting results for all RLQs, FSRQs, and SSRQs are 
listed in Table 4, where the bottom row is the Ly—Lyy rela- 
tion for the comparison RQQ sample from Lusso & Risaliti 
(2016). We first check whether the fitting results for all RLQs 
are consistent with the hypotheses of the models. Model I 
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Figure 11. Left: A *from above" view of Fig. 10, where the plane is filled by colors that represent the difference between best-fit Model I 
and Model II. The distribution of the optically selected RLQs on the LsGuz—Ly5994 plane with overlaid constant-log R lines (dash-dotted) 
is also shown. The dashed line (orange) is orthogonal to the constant-log R lines and passes through the median values of L,,4,4 and 
Lsagz of the sample. The profile (side view) of Model II at the position of the dashed line is shown in the right panel. Right: A “side 
view” of Fig. 10, where the x-axis is log R, and the y-axis is log Laxey normalized to the prediction of Model I. Therefore, Model I is the 
horizontal line at zero. The dashed curve (orange) is the profile of Model II, the top view of which is shown in the left panel. 'The blue 
dots and downward arrows are the detected and non-detected RLQs in X-rays. We calculate the running medians (purple) to smooth 
the raw data points. The running medians are apparently more consistent with Model I than Model II. 


results in Tyy = Yuv + an = ie which is 1.60 


(p = 0.05) away from y = 0.63+0-02, Due to the fact that the 
resulting error bar is relatively large for Model II, compar- 
ing Iw = Yuv = 0.72+0-08 and y results in a similar statistical 
significance (i.e., 1.50; p = 0.06). The normalization factor 
A = 0.444105 is 2.20 (p = 0.02) away from the 0.33001 
value applicable for RQQs. However, these two parameters 
together indicate that the corona component of RLQs is in- 
consistent with that of RQQs at a3.1o (p = 9x 107+) signifi- 
cance level. Consequently, RLQs cannot be treated as radio- 
quiet central engines plus jet cores in X-rays.?! The value of 
Tov = Yuv + pin dio = 0.66+0-02 resulting from Model III is con- 
sistent with y. Therefore, the physical picture of Model II 
is not supported, while the results of Model III are favored. 
We note that the only difference between the fitting results 
of Model I and Model III is that in Model III there is a 
minute amount of jet-linked X-ray emission (B = 0.03+9-01), 
This component only deviates from zero at a 1.90 (p = 0.03) 
significance level and makes up B/(A + B) S; 3% of the total 
X-ray luminosity. Note that we can use B/(A + B) to as- 
sess the typical jet contribution because the luminosities are 
normalized to values that are near to their medians before 
model fitting. We thus omit the difference between Model I 
and Model III and focus on comparing them with Model II. 


Secondly, the corresponding information criteria are 


20 To take account of the covariance of the errors of Yuy and Y radio: 
Markov chains resulting from model fitting are used to estimate 
the uncertainty of Ty, which is not necessarily larger than those 
of yay and y; uo: 
?! [n fact, the fitting results of Model II indicate that the coronae 
of RLQs have elevated X-ray activity compared with those of 
RQQs at given disk optical/UV luminosity, which motivates the 


proposal of Model III. 


listed in Table 5. For the results for all RLQs, Model III has 
smaller ICs and is preferred over Model II with AAIC = -7.61 
and ABIC = -7.61. Note that we performed another fitting 
test in which the corona term of Model II is fixed to the 
X-ray luminosity of RQQs, explicitly. The ICs of Model II 
become worse, and Model II is defeated by Model III with 
AAIC = -22.38 and ABIC = -13.19. Therefore, it is highly 
unlikely that the corona component of RLQs is identical to 
that of RQQs. Model II is as strongly disfavored if we com- 
pare it with Model I as well. 

We plot in Fig. 10 the best-fit Model I (blue), Model II 
(orange), and observational data for RLQs in the Loyey- 
L599 4 -L5anz space, in order to provide an overall geometri- 
cal impression of the best-fit models. Model III is not shown 
as its fitting results are similar to those of Model I. With axes 
in logarithmic scales, Model I is a plane while the surface of 
Model II shows upward curvature. In Fig. 11 (left) we show 
a “from above" view (i.e. looking down along the log Loyey 
axis) of Fig. 10, where the L;sooA LsGnuz plane is filled by 
colors that represent the difference between the two models. 
The distribution of RLQs overlaid with constant-R lines are 
also shown. It is easily noticed that there is a rough con- 
sistency between the direction of the color boundaries and 
the constant-R lines. In other words, the difference between 
Model I and Model II has a strong dependence on R. 

Fig. 11 (left) cannot tell us which model is better. How- 
ever, it hints at an appropriate azimuth angle for a use- 
ful “side view” that can—we should look in the direction of 
constant R. Such a side view is plotted in the right panel 


22 Model II and Model III (with fixed Yradio) have the same num- 
ber of free parameters. In this case, AAIC = ABIC and they indicate 
that Model III has a larger Lmax than that of Model II. 
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of Fig. 11, where the y-axis is normalized to the predic- 
tion of Model I. Therefore, Model I is a horizontal line at 
zero. We choose a position that is indicated by the dashed 
line (orange) in the left panel of Fig. 11 for the profile of 
Model II. This dashed line is perpendicular to constant-log R 
lines and passes through the median values of log Lasoo and 
log Lsggz. This profile of Model II is plotted as the dashed 
curve (orange) in the right panel of Fig. 11. The small dots 
and downward arrows in Fig. 11 (right) represent detec- 
tions and non-detections in X-rays, respectively. We use a 
running-box filter with à minimum of 40 and a maximum 
of 120 data points (sliding from left to right) to smooth 
the considerable scatter of the points. The step size of the 
running-box filter is 40 data points. We calculate the median 
of the data points within the box after each step. The result- 
ing running medians with uncertainty estimates are shown 
as cross symbols with error bars (purple). The running me- 
dians are consistent with zero, supporting the superiority of 
Model I in describing the data. There are systematic discrep- 
ancies between Model II and the smoothed data (running 
median) at the smallest and largest radio-loudness parame- 
ters. The latter case has the largest discriminatory leverage 
since the running medians here have relatively small error 
bars. Therefore, the superiority of Model I is apparently re- 
vealed by the median properties of RLQs with similar radio- 
loudness.?? The dispersion (0.29 dex at a 1o level) of individ- 
ual RLQs in Fig. 11 (right) has multiple sources. We assume 
typical uncertainties of 20%/30%/40% for the radio/UV /X- 
ray luminosities (e.g. Miller et al. 2011), including obser- 
vational uncertainties and quasar variability, the latter of 
which probably dominates (e.g. Gibson et al. 2008; Gibson 
& Brandt 2012). Subtracting those effects, about 0.23 dex of 
residual dispersion might be attributed to unaccounted for 
physical processes, which is similar to the magnitude esti- 
mated for RQQs (e.g. Salvestrini et al. 2019). 


23 The same conclusion is revealed if the y-axis is instead nor- 
malized to the best-fit Model II. 
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Table 4. Model fitting results for RLQs and RQQs, both of which are optically selected. 


Sample* Model Parameters 
uv radio Vea io uv radi 
I: log Lokev = @ + Yuy log Lasoo + Yagi 108 Lecute II: log Laev = log (AL? . + BL2t III: log Laev = log (ALE L7™ + BLIati) 
af 107 Yuv Y radio A B Yuv? Yradio A B? Yuv Yradio® Yaso 
0.01 0.02 0.03 0.02 0.05 0.05 0.07 0.07 0.03 0.01 0.03 0.03 
All RLQs -0.09+9:0! / 0.81+0:02 0.47003 0.227592 0-440035 0-312005 — 0.7279 0-552006 9-762 0:04 0.0399, 0-4904 — 1.00 0.181003 
0.02 0.03 0.05 0.02 0.05 0.09 0.08 0.08 0.06 0.03 0.05 0.04 
FSRQs -0.06+9:02 / 0.88+0:03 045^ 0 04 0.2659, 0-452006 — O34" yo — 0.67795 — 0.66 s — 0.76'5 9; — 0.067 — 0.49" 05 — 1.00 0.181004 
0.02 0.03 0.04 0.03 0.06 0.04 0.05 0.02 0.06 0.01 
SSRQs* -0.14°9.0 [ 0-727089 0.397504 0.227593 0-421005 0.267 9°95 0.63 0.487504 9-722 0:04 < 0.02 0.397593 1.00 0.2159: 
log Loxey = B +y log Lasoo = log (10°27. 4] 
BJ 10° y 
: E 0.01 0.01 0.02 
RQQs 0.488-01 / 0.33+0:01 0.637505 


a Note that the luminosities of the quasars are normalized to log Laxey = 27, log Lasoo = 30.5, and log Lsguz = 33.3 before model fitting (see 8 3.2.2). Effectively, radio-loudness 
parameters are normalized to log R = 2.677. The resulting normalization factors cannot be compared directly with those of other works using different normalization. 

P For the results of RLQs using Model I and the results of RQQs, the normalization factors are given in both logarithmic (o) and physical (10%) units. 

© We have fixed yuy = 0.63 of Model II for SSRQs and fixed Yradio = 1 of Model III for all RLQ samples (i.e. the underlined cells). See Appendix A for details. 

d If B is consistent with zero, we only show the upper bound of the 9096 probability interval. 

* 'The fitting results for SSRQs using Model I are consistent with the corresponding results listed in Table 7 of Miller et al. (2011), although our constraints here are 
quantitatively tighter. Note that, however, our interpretation of the results is different from that of Miller et al. (2011). See § 3.2. 


9I 


10 19 NYZ WB 


ME All RLQs ME FSRQs 


Frequency 


20 40 60 80 20 


Fret (%) 


Figure 12. The distributions of fet for all RLQs and FSRQs. 
The height of the highest bin is fixed to be identical in the two 
panels. In both cases, fet is clustered around very small values, 
while FSRQs have a long tail that extends to fet ~ 80%. 


3.2.4 The results for FSRQs and SSRQs and the amount 
of X-ray emission from jets 
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Table 5. Information criteria of different functional models for the RLQ samples. The corresponding model fitting results are listed in 
Table 4. 


Sample Information Criteria? Significance? 
Model I (k = 3) Model II (k =4 or 3)° Model III (k=4)° AAIC ABIC 
AIC BIC AIC BIC AIC BIC 

All RLQs - - —994.56 —971.61 — 1002.17 -979.22 -7.61 -7.61 

FSRQs - - —480.78 —460.90 —485.79 —465.91 -5.01 -5.01 

SSRQs —534.62 -519.67 -533.41 —518.46 - - -1.21 -1.21 


a We have bold-faced the smallest AIC and BIC of each row. The difference between the fitting results 
of Model I and Model III is always very small for all samples. We omit Model I for all RLQs and 
FSRQs, in which a small amount of jet component is present (see 8 3.2.3 and § 3.2.4). We omit 
Model III for SSRQs where the jet component seems absent (see 8 3.2.4). 

P AAIC and ABIC are calculated as the smallest values of AIC and BIC subtracted by the second 
smallest ones. 

* 'The number of free parameter is reduced by one if a parameter of Model II or Model III is fixed. See 
Appendix A for details. 
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Figure 13. Left: Fitting the data for SSRQs using Ly « Radio. The dots and downward arrows represent the detected and non-detected 
objects in X-rays. The color of each symbol represents the optical/UV luminosity, as indicated by the color bar on the right-hand side. 
We also divide the data into two bins of equal size at the median optical/UV luminosity (log L459 4 730-35) and fit them separately. 


Right: The same as the left panel for Ly œ L "to 


radio ' 


Table 6. Mean and percentiles of fet (in units of percent) for 
different RLQ samples. fjet represents the fraction of the X-ray 
emission from jets. 


Sample Mean (96) Percentile (96) 
10th 50th (Median) 90th 


All RLQs 5.2 0.3 2.3 14.4 
FSRQs 9.5 0.6 4.2 27.5 
SSRQs? - - - 5.9 


a We only list the 90th percentile for SSRQs, the X-ray lumi- 
nosities of which contain a negligible jet component. 
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Figure 14. The dependence of A on log R for RQQs and RLQs. A is the normalization factor of the coronal component (Eq. 7). A vertical line (dotted black) has been used to separate 
RLQs and RQQs. The small dots and downward arrows in blue are X-ray detected and non-detected SSRQs, respectively. The generalized Spearman rank correlation coefficient (o) 
and corresponding p-value are given in the lower-right corner. These SSRQs are further grouped into three R bins of comparable sizes, where medians with uncertainty estimates are 
calculated (large blue dots with error bars). The RQQs from Laor & Behar (2008) are small squares and leftward arrows in black, while binned data for those RQQs are large squares. 
Argo is represented by the horizontal (dashed black) line at unity, as no strong dependence on log R is reported in the literature or supported by the RQQ sample of Laor & Behar 
(2008). Note that the solid blue line is derived from the fitting results for SSRQs using Model I (see Table 4). The shaded regions represent lo: uncertainties, which are small for the 
case of RQQs. Note that the good match between dashed and solid lines at log R = 1 is not an enforced prerequisite but a direct result of the model fitting. 
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FSRQs and SSRQs are RLQs that are generally ob- 
served at relatively small and large inclination angles, re- 
spectively (e.g. Urry & Padovani 1995). If part of the X-ray 
luminosity is from the beamed core region of jets, the contri- 
bution of this component is expected to be larger in FSRQs 
than SSRQs. Indeed, their properties seem to be different as 
suggested in 8 3.1; FSRQs are generally more X-ray lumi- 
nous than SSRQs at fixed Luv and R (see Fig. 6). 

For the fitting results for FSRQs in Table 4, the pa- 
rameters of Model I do not support Tuy = y, while those 
of Model II do not support a component that is consistent 
with the coronae of RQQs. Model III might point to the 
most plausible scenario with Tuy = yuv + Vea dio = 0.67+9-04, 
Note that the jet component is set to zero in Model I, while 
its fraction in Model III is still only B/(A+ B) « 7.0%. There- 
fore, we do not distinguish between those two models, similar 
to the case of 8 3.2.3. Model II is disfavored in comparison 
with Model III (Model I) with AAIC = —5.01 (-2.27) and 
ABIC = -5.01 (-6.24). Note that the fitting results for FS- 
RQs are largely consistent with those for the sample of all 
RLQs in $ 3.2.3, except that the Bp term is more sig- 
nificant. Note that the residual scatter after subtracting the 
effects of measurement errors and variability is about 0.17 
dex for FSRQs (and SSRQs below). 

As for SSRQs, the fitting results of Model III are almost 
identical to those of Model I with B consistent with zero (see 
Table 4). In this case, B is totally redundant in that it in- 
creases the number of free parameters without improving the 
performance. We thus omit Model III for SSRQs. Model I 
results in Tuy = Yuv + Yiadio = 0.619-03, consistent with y. 
The fitting results of Model II do not support a corona com- 


ponent that is identical to those of RQQs (A = 0.42+0-06 vs. 


108 = 0.33+9-01) and Yradio = 0.48+0.06 departs strongly from 
unity. Furthermore, Model II is not favored in comparison 
with Model I with AAIC = -1.21 and ABIC = -1.21 (see Ta- 
ble 5). Note that the model selection results here are not as 
significant as for the cases of all RLQs and FSRQs, due to 
a smaller sample size and lower X-ray detection fraction for 
SSRQs. 

Benefiting from its flexibility, Model III either is se- 
lected as the best model or can be consistent with the best 
model for all samples. Using the results of Model III in Ta- 
ble 4, we now assess the importance of the jets in account- 
ing for the X-ray luminosities of RLQs. We input Loi 
and Lsgyz as well as estimates of parameters into Model III. 
We quantify the fraction of jet emission for each RLQ using 


JAN Yradio Yiadio Yuv Yradio 
fiet = BL GHz (ALt Esco +B LsGHy) The mean and three 


percentiles (10, 50, and 90th) of fjet for samples that might 
have a jet component are listed in Table 6, while the full 
distributions are plotted in Fig. 12. For SSRQs, we only list 
their 90th percentile, which indicates that for 90% of SSRQs 
the jet component contributes < 6% of the observed nuclear 
X-ray emission. Note that B/(A- B) is generally between the 
mean and median of fjet. 

The distributions of fjet are highly asymmetric in 
Fig. 12. As a consequence, the median is always smaller than 
the mean in Table 6. The median is probably a better statis- 
tic to assess the typical value of fjet, which is < 5% for all 
samples. FSRQs are the group of RLQs that have the largest 
jet contribution. Notably, the 10th and 90th percentiles span 
several decades in Table 6, which indicates a large variance 


of fjet from object to object. The mean is strongly affected 
by a small number of RLQs that have relatively large fjet. 
Less than 10% of FSRQs have fjet Z 30%, which manifest 
themselves as a long tail in the right panel of Fig. 12. Thus, 
only in a minority of FSRQs does the jet component become 
important. 

At fixed Lasoo A: RLQs are generally a factor of 1.5-4 
times more X-ray luminous than RQQs, depending on the 
radio-loudness parameter (see Fig. 6). Given our limits on 
fet, this difference is mainly caused by the difference between 
the coronae of RLQs and those of RQQs, not the emission 
from jets. Furthermore, typical FSRQs are more X-ray lumi- 
nous than SSRQs at fixed Ly, and R by 1096-2096 (Fig. 6), 
which is only partially attributed to the jet component ac- 
cording to our estimation here, given that the median fjet is 
only 4.2% for FSRQs. Therefore, other processes are proba- 
bly involved as well (e.g. a small level of jet contribution to 
the optical/UV emission for FSRQs) to explain the differ- 
ence between FSRQs and SSRQs. 


3.3 The corona-jet connection in SSRQs 


Owing to the apparent absence of significant X-ray emis- 
sion from jets in general, SSRQs form a useful sample to 
investigate the corona-jet connection as implied by Model I. 


In § 3.2.1, we interpret Model I as Ly « RY catio 12°" radio 


instead of Lx o res LY, The argument there was based on 
the fact that yuv *Ya dio 18 close to y (the slope of the Lx—Luy 
relation for RQQs), which might suggest an identical relation 
between ALx/ Ly and ALyy/Lyy for both the RQQ and RLQ 
populations. The results of § 3.2.3 and 8 3.2.4 indeed con- 
firm this approach. We now show that the former expression 
is preferred even if the properties of RQQs are not consid- 
ered, and our interpretation of the corona-jet connection is 
thus independent of outside assumptions. We fit the data 


Vaas 
5GHz’ 
= 0.36+0-02 respectively. 


for SSRQs using Loyey « Rio and Lopey « L which re- 


sults in Yl dio = 0.22*0-04 and y’ 


.—0.04 © radio. —0.02* 
The former value is consistent with that from the three- 
parameter model in Table 4 (i.e. "adio = 0.320-03 in Model I 


for SSRQs). The results are shown in Fig. 13, where the 
data points are color-coded by their optical/UV luminosi- 
ties. In the left panel of Fig. 13, the gradient of data-point 
colors appears perpendicular to the line of Loyey « RY ratio, 
while in the right panel Loyey œ joe is not able to sep- 
arate high-luminosity and low-luminosity objects. To con- 
firm this, we divide SSRQs into high-luminosity and low- 
luminosity bins at the median log L sọọ = 30.35 and fit 
them separately. The results are shown in Fig. 13 as well. 
The lines in the left panel are largely parallel, with the high- 
luminosity and low-luminosity fits well straddling the orig- 
inal fit. The right panel is complex, and no useful informa- 
tion can be easily extracted. Furthermore, reading Fig. 13 
(left) and Fig. 8 (bottom) together reveals an attractive 


+ 
Yuv VY radi o 


property: Ly « RY radio LP can be decomposed into two 


statistically independent relations such that Ly œ Radio and 


Ly « AM ^ Therefore, Lx depends primarily on Luy and 


R, while its relation with Lradio is a by-product. Miller et al. 


24 The commonly assumed situation where the properties of RLQ 
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(2011) notice that Lx/Lx goo has a stronger dependence on R 
than Lradio (especially for SSRQs), which is consistent with 
the idea that R is the more meaningful parameter for de- 
scribing the corona-jet connection of RLQs. Consequently, 
the index of the X-ray-optical/UV relation for the corona is 
universal among RLQs and RQQs (i.e. Tuv = yq, tYuv = Y); 
which is now a direct result of separated correlation analysis 
of two quasar populations. 

More importantly, due to the above exposed statistical 
independence, another process is probably in control of the 
corona-jet connection in RLQs aside from the disk-corona 
interaction that is at work for both RLQs and RQQs. Follow- 
ing § 3.2.1, the discussion is continued in the context of Eq. 7 
where A is a function of R. We plot A of SSRQs divided by 
that of RQQs against R in Fig. 14. The scattered individual 
data points together with binned data points (separated by 
the 33rd and 66th log R percentiles of SSRQs) are shown. We 
calculate the generalized Spearman’s rank correlation coeffi- 
cient of the individual points in Fig. 14, using the Astronomy 
Survival Analysis package (Lavalley et al. 1992). The corre- 
lation is strongly supported with a p-value of x 1x 1079. The 
solid blue line (with shaded region showing estimated uncer- 
tainties) in Fig. 14 is a projection of Model I in Table 4 for 
SSRQs, and not a new fitting of the individual data points. 
The normalization factor of RQQs (i.e. Argo) is at unity 
(dashed black line), where the shaded region represents its 
small uncertainty. We are unaware of any established radio 
dependence of the Ly-L,y relation for RQQs (e.g. Steffen 
et al. 2006; Just et al. 2007; Lusso et al. 2010). More im- 
portantly, the radio luminosity of RQQs might not be dom- 
inated by the synchrotron emission from quasar jets (e.g. 
Kellermann et al. 2016; Panessa et al. 2019). After subtract- 
ing the effects of measurement errors and quasar variability, 
the physical intrinsic scatter of the Lx—Lyy relation is about 
0.2 dex (e.g. Lusso & Risaliti 2016; Chiaraluce et al. 2018; 
Salvestrini et al. 2019). There is thus limited error budget 
for additional physical processes involving radio luminosity 
(or radio-loudness parameter) to play an important role in 
RQQs as in RLQs. We show the 59 RQQs from Laor & Be- 
har (2008) as black symbols, where they are also grouped 
into two R bins. These PG quasars are consistent with the 
horizontal line and do not show an R dependence. 

'The dashed and solid lines are satisfyingly consistent 
with each other at R « 10; this agreement has not been 
enforced but arises entirely from the independent model fits. 
Consequently, if we express the X-ray luminosity of RQQs 
as 


Lx = Arggliv: (10) 


the X-ray luminosity of the corona of RLQs can be written 
as 


Lx = Aggg(R/10)" 4° LT. . (11) 


The good match between Argo and Arro around R = 10 
(instead of some much smaller radio-loudness parameter) in 
Fig. 14 prevents extension of the corona-jet connection of 
RLQs into the radio-quiet regime. 

It is also tempting to attempt to connect RQQs and 


coronae are identical to those of RQQ coronae probably needs 
finely tuned jets to produce such a property. 
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RLQs with a single function that contains a “jet component”, 
which grows smoothly from being negligible at R < 10 to 
being significant at R > 10. However, this hypothetical com- 
ponent unavoidably has a strong dependence on optical/UV 
emission (e.g. Fig. 8) and a relatively weak dependence on 
radio emission (e.g. Fig. 13), in contrast with expectations 
for a “jet component”. We thus suggest that a break at R ~ 10 
is probably required, which points to an intrinsic difference 
between the innermost accretion flows of RLQs and RQQs. 
Therefore, the radio-loudness parameter of quasars is more 
fundamental than simply empirical. Future works focusing 
on the X-ray properties of quasars with 0 < log R < 2 might 
reveal the exact transition clearly. 


4 DISCUSSION 
4.1 Comparison with literature results 


Given that we find the corona is still responsible for most 
of the X-ray emission of RLQs (as for RQQs), it is of value 
to compare with literature results on the spectral, imaging, 
and timing properties of radio-loud AGNs in X-rays to assess 
overall consistency. 


4.1.1 X-ray spectral properties 


X-ray spectral analyses of radio-loud AGNs, except for some 
FSRQs, using high-quality X-ray data rarely reveal a distinct 
jet component. For example, the X-ray spectra of BLRGs 
like 3C 382 (e.g. Gliozzi et al. 2007; Ballantyne et al. 2014; 
Ursini et al. 2018), PKS 2251+11 (e.g. Ronchini et al. 2019), 
3C 120 (e.g. Rani & Stalin 2018), and 3C 390.3 (e.g. Sam- 
bruna et al. 2009; Lohfink et al. 2015) are consistent with be- 
ing disk/corona-related. Specifically, their observed primary 
power-law, reflected Fe Ka line (and Compton hump), and 
soft excess are all established features common to radio-quiet 
AGNS as well. Also, importantly, those studies utilizing NuS- 
TAR data in the hard X-rays generally find a high-energy 
cutoff at ~ 100 keV (e.g. Ballantyne et al. 2014; Lohfink et al. 
2015; Rani & Stalin 2018), signifying a thermal Comptoniza- 
tion process of X-ray continuum production as for radio- 
quiet Seyfert galaxies. Lohfink et al. (2017) analyzed the 
Swift/ NuSTAR X-ray spectra of an SSRQ (4C 74.26); no 
sign of jet emission but a high-energy cutoff of 183+31 keV 
was revealed. 

Note that a jet-dominated spectrum usually shows an 
“upturn” (in SED representation) in hard X-rays, extend- 
ing to MeV energies (and sometimes beyond), instead of 
a “rollover”. Indeed, the upturn is present in those FSRQs 
where the jet emission dominates (e.g. Paliya et al. 2016; 
Ghisellini et al. 2019) or has a significant contribution (e.g. 
Grandi & Palumbo 2004; Madsen et al. 2015) in X-rays. 
However, those objects with strongly beamed jet X-ray emis- 
sion are only a small portion of their parent population (see 
§ 2.1.4). Some sample-based archival X-ray spectral studies 
finding flat X-ray power-law continua for RLQs relative to 


25 The reflection features of BLRGs in statistical samples are 
known to be weaker than those of radio-quiet Seyfert galaxies, 
which can be attributed to mechanisms other than the dilution 
by continuum X-ray emission from jets (e.g. Eracleous et al. 2000). 
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Figure 15. The distributions of RQQs (open squares) and SSRQs (solid circles) in the L599 ~@ox plane. The RQQs are from Laor & 
Behar (2008). SSRQs are color-coded according to their radio-loudness parameters (following the color-bar on the right-hand side). The 
dashed (black) line is the a@x—Luy relation for RQQs, derived from the Ly-Lyy relation in Fig. 4. The solid (light blue) line is derived 
from the fitting results for SSRQs using Model I (see Table 4), fixed at the median radio-loudness (log R = 2.69). We interpret this plot 
as a HID for quasars, in analogy to that for BHXRBs. We propose that the @ox—Lyy relation for RQQs corresponds to an approximate 
“jet line" for quasars; the quasars lying well to the right of the jet line have powerful relativistic jets, which are quenched for quasars 


that are “on” or to the left of the jet line. 


those of RQQs (e.g. Wilkes & Elvis 1987; Lawson & Turner 
1997; Reeves et al. 1997; Page et al. 2005) utilize heteroge- 
neous RLQ samples, dominated (at the 60-100% level) by 
flat-spectrum RLQs, and potentially extreme objects. Thus, 
these samples cannot be reliably used to constrain the X-ray 
spectral properties of general RLQs, the focus of this work. 
Generally, once FSRQs are excluded, the X-ray spectra of 
the remaining BLRGs, narrow-line radio galaxies, and SS- 
RQs are similar to those of radio-quiet AGNs (e.g. Lawson 
et al. 1992; Galbiati et al. 2005; Grandi et al. 2006).26 
Gupta et al. (2018) investigated luminous radio galaxies 
and their radio-quiet counterparts with comparable Edding- 
ton ratios and black-hole masses. The radio galaxies were 
found to be, on average, a factor of ~ 2 more luminous in 
hard X-rays than the radio-quiet AGNs, while their power- 
law spectral slopes and high-energy breaks were similar. The 
authors concluded that the X-ray emission in both samples 
is produced by à common mechanism. However, the X-ray 
production efficiency is apparently higher for the radio-loud 


26 Note that even if the X-ray spectra of general RLQs are some- 
what flatter than those of RQQs, this would not necessarily 
demonstrate that the X-ray emission is jet linked (e.g. Laor et al. 
1997). 


group, which is consistent with our suggestion of a corona- 
jet connection (e.g. 8 3.3). Gupta et al. (2020) extended the 
comparison and showed that Type 1 AGNs are not X-ray 
louder than Type 2 AGNS, within both radio-loud and radio- 
quiet groups. Therefore, the dependence of observed X-ray 
luminosity on viewing angle seems weak, inconsistent with 
beamed jet X-ray emission. 


4.1.2 X-ray imaging properties 


The Chandra observatory with sub-arcsec angular resolu- 
tion has detected kpc-scale X-ray jets from many RLQs (e.g. 
Harris & Krawczynski 2006), which might contribute to the 
X-ray fluxes we use (see 8 2.1.2). The observed X-ray lu- 
minosities of these extended jets are generally only a few 
percent those of the cores (e.g. Marshall et al. 2005, 2018), 
consistent with the amount of jet contribution we estimated 
in 8 3.2.4. 

Imaging the nuclear region of quasars in X-rays is cur- 
rently limited to indirect methods (e.g. gravitational lens- 
ing). The sizes of the nuclear X-ray emission regions derived 
from gravitational-lensing studies seem to be systematically 
larger in RLQs than in RQQs, for a given SMBH mass, al- 
though the source statistics are presently very limited (e.g. 
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Figure 16. The SSRQs of this paper and RQQs of Laor & Behar (2008) in the Lyagio-Luy plane. The symbol colors represent observed 
X-ray luminosities divided by those expected for RQQs at given matched optical/UV luminosity, Lx/Lx,RQQ, as indicated by the color-bar 
on the right-hand side. X-ray non-detected SSRQs (11.3%) are not shown, the distribution of which is consistent with the detections (see 
Fig. 14). Radio non-detected RQQs are shown as downward arrows. Four lines indicate log R = 1 to 4, from light to dark. There is an 


apparent correlation between Lx/L. goo and R for SSRQs (cf. Fig. 6). 


Burak Dogruel et al. 2019). Therefore, we might not expect 
the coronae of RLQs and RQQs to have the same physical 
properties, the former of which could be larger in units of 
black-hole radius. 


4.1.8 X-ray timing properties 


Long-term timing studies are generally more observation- 
ally intensive than single-epoch spectral studies of AGNs in 
X-rays, and thus are sparse for general radio-loud AGNs.77 
Leighly et al. (1997) investigated the 9-month X-ray variabil- 
ity of the BLRG 3C 390.3 and found its fractional amplitude 
of variability is about 33%. This amount of variability is not 
apparently different from those of radio-quiet Seyfert galax- 
ies given the luminosity of 3C 390.3 and the timescales this 
study covers. Furthermore, the hour-to-year power spectra of 
several BLRGs (e.g. 3C 120, Marshall et al. 2009; 3C 390.3, 
Gliozzi et al. 2009; 3C 111, Chatterjee et al. 2011) are also 
similar to those of radio-quiet Seyfert galaxies. The X-ray 


27 Those rare radio-loud AGNs with one of their jets pointing 
close to our line of sight have generally gained more attention 
(e.g. Marscher et al. 2010; Hayashida et al. 2015). However, their 
X-ray variability probably has a different origin and is strongly 
affected by Doppler boosting effects. 
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variability properties of RLQs are poorly constrained for sys- 
tematically derived samples as well, let alone for subsamples 
with established radio slopes. Sambruna (1997) investigated 
the soft X-ray variability of a few FSRQs on month-to-year 
timescales (> 6 epochs) and found a typical variability am- 
plitude of 10-3096, roughly consistent with that of RQQs 
(e.g. Gibson & Brandt 2012). Gibson & Brandt (2012) stud- 
ied the X-ray variability of about 20 non-BAL RLQs and 
found suggestive evidence that they are less variable than 
RQQs. However, the observations of these RLQs generally 
have only 2-3 epochs. Future X-ray variability studies utiliz- 
ing RLQ samples with measured radio spectral slopes and 
observations with more epochs might improve our under- 
standing of the X-ray variability properties of RLQs and 
shed further light on the origin of their X-ray emission. 


4.2 A jet line for AGNs 


In this section, we investigate further the connection be- 
tween the disk/corona and jets from the perspective of the 
Qox—Luy relation. Fig. 15 shows the transposed aox-Luy plot 
for RQQs and SSRQs, which is no more than an alternate 
representation of Fig. 8 (bottom). However, we gain new 
insights from it by analogy with BHXRBs. 

We interpret Fig. 15 as a quasar-version of the hardness- 
intensity diagram (HID). The HID is frequently used in 
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BHXRB studies. The y-axis serves as a proxy for the ac- 
cretion rate of the black hole, while the x-axis represents 
the relative contributions of the thermal (at 2500 À) and 
power-law (at 2 keV) components that are radiated from 
the accretion disk and corona, respectively. As discussed in 
detail earlier, a significant jet contribution in X-rays for SS- 
RQs is highly unlikely, legitimizing our approach. Note that 
the data points in an HID for BHXRBs often represent the 
evolutionary stages of a single black hole, which in Fig. 15 are 
represented by an ensemble of black holes that span ranges 
in accretion rate, ^hardness ratio", and jet power. 

In Fig. 15, the dashed line (black) represents the @ox—Luy 
relation for RQQs, which was transformed from the Lx—Lyy 
relation in Table 4 (bottom row) and Fig. 4. For SSRQs, we 
utilize the fitting results of Model I in Table 4 to derive their 
Q@ox—Luy relation (light blue) at fixed log R = 2.69, which is 
the median radio-loudness parameter for SSRQs in our sam- 
ple. The dashed and solid lines are parallel to each other, as 
a consequence of luy « y. Open squares and filled circles are 
RQQs from Laor & Behar (2008) and SSRQs, respectively. 
SSRQs are further color-coded by their radio-loudness pa- 
rameters. Among RQQs, more luminous objects have steeper 
Qox, and hence a smaller power-law fraction. SSRQs to first 
approximation follow this trend as well. However, the radio- 
loudness parameter also plays an important role such that 
those more radio-loud quasars have flatter vo, at given Luv, 
which is a consequence of Agta « ZigoQ(R/ 10) "radio (as per 
Eq. 10 and Eq. 11). To summarize, as shown in Fig. 15: 


(i) The most radio-loud SSRQs (say, with log R = 4) lie on 
a line near the right edge of the data points that is parallel 
to the a@x—Luy relation for RQQs. 

(ii) Horizontally shifting such a line leftward, it passes 
through SSRQs with generally decreasing R. 

(iii) Finally, the line arrives at the region populated by 


RQQs. 


Radio observations of BHXRBs suggest the existence of a 
critical line in the HID such that jets are active on the 
right-hand side of this line, but seem to be quenched once 
a BHXRB moves leftward and passes the line (e.g. Fender 
et al. 2004). Fig. 15 might reveal a similar jet line for AGNs, 
which is approximately the a@ox—Lyy relation for RQQs. The 
jet line for BHXRBs seems to be tilted in a way such that 
the quenching happens at a larger hardness ratio when the 
luminosity is lower (e.g. Fender et al. 2009), similar to the 
alignment of the black-dashed line in Fig. 15. 

The jet line adds to the mounting evidence for the 
idea that black-hole accretion/jet physics is largely scale- 
invariant (e.g. Merloni et al. 2003; McHardy et al. 2006; 
Arcodia et al. 2020). RLQs and RQQs are likely scaled-up 
versions of BHXRBs in different accretion states (e.g. Nipoti 
et al. 2005; Kórding et al. 2006; Sobolewska et al. 2009). 


4.3 Origin of radio-loudness 
4.9.1 The overall corona-disk-jet relation 


We show in Fig. 16 the SSRQs in the Lyagip-Luy plane, re- 
orienting the discussion from, hitherto, the nature of X-ray 
emission to the origin of powerful radio jets. The symbol 
colors represent the X-ray luminosities divided by those of 
RQQs at given matched Luy, i.e. Lx/Lx RQQ, utilizing the 


Lx-Lyy relation in Fig. 4. The X-ray non-detections (11.3%) 
are omitted in Fig. 16, the distribution of which is similar 
to that of the detections for SSRQs (cf. Fig. 14); omitting 
non-detections will not bias the result. T'he 59 RQQs of Laor 
& Behar (2008) are also shown for comparison, where radio 
non-detected objects are represented by downward arrows. 
Apparently, RQQs follow a different track in Fig. 16 than 
SSRQs such that typical RQQs are about 3 orders of mag- 
nitude less radio-luminous than RLQs at a given accretion 
rate (e.g. Rawlings 1994). Whether the dichotomy between 
RLQs and RQQs is real or not has long been debated (e.g. 
Kellermann et al. 1989, 2016; Ivezić et al. 2002; Cirasuolo 
et al. 2003; Balokovié et al. 2012; Condon et al. 2013). The 
results of this paper (e.g. Fig. 14) support the idea that 
RLQs and RQQs are intrinsically different objects, regard- 
less of the bimodality of the distribution of R (e.g. Padovani 
2016). The accretion flow for RQQs is probably in a state 
that is not compatible with launching powerful relativistic 
jets. 

The most prominent relation for SSRQs in Fig. 16 is 
that between Lradio and Luv, which is not one of the main 
topics of this paper but reveals probably no less important 
physics (e.g. Serjeant et al. 1998; van Velzen & Falcke 2013). 
The Lradio—Luv relation is expected from the near-linear cor- 
relation between the radiative powers of the outflowing jets 
and the infalling accretion disk for RLQs (e.g. Ghisellini 
et al. 2014), supporting a dependence of the jet power on 
accretion rate (e.g. Rawlings & Saunders 1991). Therefore, 
in addition to the corona-jet connection investigated in this 
work and disk-corona interplay, RLQs also feature a disk- 
jet connection. For the first time, not only has the overall 
corona-disk-jet relationship of RLQs been revealed, but also 
the number of responsible physical processes has been deter- 
mined. In short, the corona-disk-jet relationship is three-fold 
such that any two of these structures are physically linked.?8 
Furthermore, the underlying drivers for these three relations 
are independent, and none of these three relations is simply 
a side effect of the other two. For example, we found that 
R is a more useful property to consider than Lradio when 
the corona-jet connection is examined in 8 3.3, which now 
has a physical reason. Specifically, the radio luminosity is af- 
fected by both the disk-jet and corona-jet connections; using 
R minimizes the intrusion of the former into studies of the 
latter. Fig. 16 demonstrates more clearly the independence 
between the disk-jet and corona-jet connections.?? There is 
an apparent correspondence between R and Ly/L. goo (cf. 
Fig. 6) such that the symbol-color gradient appears largely 
perpendicular to the constant-log R lines. From Fig. 16, it 
is strongly suggested that the spread of R within the RLQ 
group is directly caused by small-scale processes in the vicin- 
ity of SMBHs (relating black-hole spin, Eddington ratio, and 
magnetic field) that, at the same time, affect the corona. 


4.3.2. The role of black-hole spin 


Models for the launching of quasar jets need to explain our 
corona-jet connection in RLQs as well. The role of rapidly 


28 RQQs share only a common disk-corona interplay with RLQs. 
7? We have already shown the independence between the corona- 
jet connection and disk-corona interplay in 8 3.3. 
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Figure 17. The cumulative distribution functions of Lx/ Ly &QQ 
for RQQs (from Lusso & Risaliti 2016; black), SSRQs (blue), and 
a subset of SSRQs with R > 1000 (red). Almost no RQQs are as 
X-ray luminous as the top 15% of all SSRQs or the top 30% of 
SSRQs with R > 1000. Similarly, almost no SSRQs are less X-ray 
luminous than the bottom 15% of RQQs. 


spinning black holes can thereby be constrained, which are 
often taken to be an important ingredient for jet launching 
(e.g. Blandford et al. 2019). At first sight, it seems plausi- 
ble that since higher prograde black-hole spins increase both 
the jet-production efficiency and radiative efficiency of the 
inner accretion flow, a correlation between R and Lx/ Lx REQ 
may be naturally produced (cf. Fig. 6 and Fig. 16). How- 
ever, the majority of spin parameters resulting from X-ray 
Fe Ka modeling of SMBHs are = 0.9 (e.g. Reynolds 2019, 
cf. Laor 2019), even for radio-quiet Seyfert galaxies.) Un- 
less some mechanism can largely decouple black-hole spin 
and the corona for (only) radio-quiet AGNs, those RQQs 
with the highest black-hole spins are then expected to be 
as X-ray luminous as their most radio-loud counterparts. In 
Fig. 17, we show the cumulative distribution functions of 
Lx/Lx,RQQ for our SSRQs and the RQQ sample of Lusso & 
Risaliti (2016). Apparently, there are almost no correspond- 
ing RQQs to the top 15% of SSRQs (i.e. the most radio-loud 
group) in Lx/Lx goQ. If we instead use a subset of SSRQs 
with R » 1000 (red curve in Fig. 17), the value can be re- 
laxed to the top 3096. 

Considering a worst-case scenario where the X-ray 
Fe Ka method significantly overestimates the spin parame- 
ter of radio-quiet AGNs and all RQQs have slowly rotat- 
ing black holes, a group of X-ray under luminous RLQs 
is expected since black holes with high retrograde spin 
parameters can likely launch relativistic jets as well (e.g. 
Tchekhovskoy & McKinney 2012), the radiative efficiency 
of which will be lower than those of RQQs. However, no 
SSRQ extends beyond the 15% of RQQs at the lower end in 
Lx[Lx RQQ (see Fig. 17). 

A decisive role of black-hole spin is also inconsistent 
with the jet line for AGNs in § 4.2. Indeed, a BHXRB can be 
in both jet-active and jet-inactive states without a significant 


9? Using the Soltan (1982) argument, constraints on the accretion 
efficiency may also imply that most RQQs have rapidly spinning 
SMBHs (e.g. Elvis et al. 2002; Yu & Tremaine 2002; Shankar 
et al. 2020). However, these Soltan-argument constraints still have 
a sizeable variance owing to the uncertainties of several utilized 
parameters (e.g. Brandt & Alexander 2015; Comastri et al. 2015). 
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change of the black-hole spin, a conclusion which by analogy 
probably applies to SMBHs as well. 

Therefore, the corona-jet connection supports the idea 
that rapid BH spin is a necessary but not sufficient condition 
for the development of jets, which likely requires a second, 
independent factor (e.g. Blandford et al. 2019). 


4.3.8 The role of magnetic flux / topology 


We suggest that the magnetic flux threading the SMBH is 
likely the key physical factor controlling the radio-loudness 
of AGNs (e.g. Sikora & Begelman 2013). Merloni & Fabian 
(2001, 2002) discussed the possibility of a connection be- 
tween the coronae and relativistic jets in the most direct 
form such that the former contains the magnetic fields re- 
quired for the launching of the latter. Therefore, only when 
these magnetically dominated coronae are sufficiently strong 
can jets be launched from them, and a possible corona-jet 
connection is expected. Indeed, our results support the idea 
that the activity of the corona is at least a tracer for the mag- 
netic flux responsible for the launching of relativistic jets of 
RLQs. Note that the topology of the magnetic field might 
play a role in determining accretion states and jet launching 
as well (e.g. Livio et al. 2003; Dexter et al. 2014). 

Furthermore, the inner accretion flows of RLQs may be 
a magnetically arrested disk (MAD; e.g. Igumenshchev et al. 
2003; Narayan et al. 2003), in contrast with the standard and 
normal evolution (SANE; e.g. Narayan et al. 2012) scenario 
for RQQs. In the MAD state, the magnetic flux thread- 
ing the black hole is regulated by the accretion rate (e.g. 
Tchekhovskoy et al. 2012; Zamaninasab et al. 2014), perhaps 
thereby naturally producing a Lyagio—Luy relation for RLQs 
(see Fig. 16). As a consequence, to explain the corona-jet 
connection established in this work, we need another factor 
that can independently affect the magnetic flux, in addition 
to the accretion rate. We speculate that the thickness of 
the MAD is a good candidate since a thicker disk is able 
to support a stronger magnetic field (e.g. Tchekhovskoy & 
McKinney 2012). The MAD disk contains optically thick ac- 
cretion streams that are rapidly spiralling toward the black 
hole and magnetically dominated low-density voids that can 
contribute to hard X-ray emission via the inverse- Compton 
scattering process (i.e. behaving like coronae; e.g. Igumen- 
shchev 2008). The magnetospheric radius (i.e. size of the 
MAD disk)?! also positively depends on the disk thickness 
(e.g. Xie & Zdziarski 2019). Therefore, a thicker disk can 
launch more powerful jets (e.g. Tchekhovskoy & McKin- 
ney 2012) as well as produce an enlarged X-ray emitting 
corona, producing the corona-jet connection. Note that the 
disk thickness might be affected by other more fundamental 
parameters (e.g. the Eddington ratio). 


5 SUMMARY AND FUTURE PROSPECTS 
5.1 The main results of this paper 


In this paper, we construct a large uniform RLQ sample 
with high-quality multi-wavelength coverage (see § 2). We 


3! Tt is convenient to treat this radius as the inner truncation 
radius of the standard thin disk outside the MAD disk. 
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investigate the X-ray-optical/U V-radio relation and the na- 
ture of the X-ray emission of RLQs. The main results are 
summarized as follows: 


(i) The X-ray luminosities of RLQs are systematically 
larger than those of appropriately matched RQQs and de- 
pend on Luy, Lradio (or R), and radio spectral index. The 
slope of the Lx-Lyy relation for SSRQs is consistent with 
that for RQQs to within z 3% (0.61+0:04 vs. 0:6310), sup- 
porting the idea that SSRQ X-ray emission is mainly from 
the corona instead of the jets. The corresponding intercept 
for SSRQs is always larger than that for RQQs (by factors 
up to ~ 3) and increases with R, suggesting a coordination 
between the activities of the jets and corona. Typical FSRQs 
are more X-ray luminous than SSRQs for given Ly, and R, 
likely involving more physical processes. See § 3.1. 

(ii) We perform model fitting and formal model selection 
in § 3.2. The models and methods are described in § 3.2.1 and 
§ 3.2.2, respectively. In the model selection for all RLQs, a 
model with a strong jet X-ray contribution is ruled out. The 
preferred models either attribute all X-ray emission to the 
corona or contain only a small amount of jet X-ray emission. 
See § 3.2.3. 

(iii) For SSRQs and FSRQs considered as groups, a model 
with a strong jet contribution to the X-ray emission is not 
favored either. The typical jet contribution for SSRQs is 
consistent with zero, again indicating that the corona is the 
dominant X-ray emitting structure. The observed X-rays 
from FSRQs might have a contribution from the jets, which 
is, however, only significant for < 10% of FSRQs. See § 3.2.4. 

(iv) Given that the jets are apparently not contributing 
substantially to the X-ray luminosities of SSRQs, we discuss 
the corona-jet connection in detail in § 3.3. Based on pat- 
terns apparent in our data, we suggest that the disk-corona 
interplay and corona-jet connection operate in largely inde- 
pendent ways. The former process is at work for both RQQs 
and RLQs, while the later process affects only RLQs. 

(v) The literature results on the X-ray spectral, imaging, 
and timing properties of radio-loud AGNs are generally con- 
sistent with the idea that the disk/corona is responsible for 
the bulk of the X-ray emission of most RLQs. See § 4.1. 

(vi) We propose that the a@ox—Lyy relation for RQQs ap- 
proximately defines a jet line for AGNs, drawing an analogy 
with BHXRBs. See § 4.2. 

(vii) RLQs feature corona-jet, disk-corona, and disk-jet 
connections, each of which seems to be driven by a distinct 
physical process. The magnetic flux threading the SMBH 
instead of the black-hole spin is most likely controlling the 
jet-launching process of RLQs. See § 4.3. 


Perhaps the most important result of this paper is that 
the jets generally contribute much less to the X-ray emission 
than was previously thought. The corona-jet connection and 
other results are then unavoidable corollaries. Our modeling 
results do allow some modest X-ray emission to arise from 
the jets in general, as is expected on basic physical grounds 
and also sometimes observed via X-ray imaging studies (see 
§ 4.1.2). However, aside from a minority of FSRQs, this jet- 
linked X-ray emission generally appears secondary relative 
to that from the corona. 


5.2 Future work 


Many lines of study that extend our results can follow as 
future work. First, as mentioned at the end of § 3.3, the 
apparent break at R ~ 10 needs further investigation, pos- 
sibly with a large sample of 0 < logR < 2 quasars. Sec- 
ondly, X-ray spectral analyses for a sample of SSRQs along 
a sequence of increasing Lx/LxrQQ (or R) might reveal how 
the corona changes, following the increasing production ef- 
ficiency of jets. Thirdly, in addition to the viewing-angle 
based unified model, an evolutionary aspect of the AGN phe- 
nomenon is likely essential as well (e.g. Klindt et al. 2019). 
The identification of a jet line makes Fig. 15 a useful diag- 
nostic for accretion states of various AGN populations. For 
example, at least some weak-line quasars and BAL quasars 
are intrinsically X-ray weak (e.g. Leighly et al. 2007; Luo 
et al. 2014), thus perhaps representing different states than 
normal quasars. Fourthly, one still-missing piece of the jig- 
saw for the X-ray properties of RLQs is systematic sample- 
based studies of their X-ray variability. It would be valuable 
to compare RQQs and RLQs (with established a; values) 
as well as AGNs and BHXRBs in this regard. Finally, the 
total radio luminosity we use traces the time-averaged jet 
power over perhaps up to ~ 107 years, while the X-ray lu- 
minosity is a more instantaneous tracer of the current ac- 
tivity of the corona. Therefore, the corona-jet connection we 
found is probably subject to a delay and variability induced 
“smearing” already. We might anticipate a closer correlation 
through long-term radio/X-ray variability monitoring of in- 
dividual RLGs, in particular the high-frequency radio emis- 
sion that is from regions closer to (or even co-spatial with) 
the X-ray emitting corona and has shorter delays. 

'These future projects will benefit from on-going and 
scheduled multi-wavelength sky surveys, especially in the ra- 
dio bands that are witnessing a resurgence (e.g. Norris 2017). 
The eROSITA telescope (Merloni et al. 2012) and Large 
Synoptic Survey Telescope (Ivezié et al. 2019) will provide 
copious time-domain X-ray and IR/optical/UV data that 
can advance our understanding of the variability properties 
of AGNs. 
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APPENDIX A: THE FIXED PARAMETERS IN 
MODEL FITTING 


We fix Yradio = 1 for Model III in Table 4 to have a better 
constraint on the X-ray contribution from the jets, which is 
one of our scientific goals. This treatment is similar to that 
of Browne & Murphy (1987). The AD term of Model III 
represents the X-ray luminosity from the core region of the 
quasar jets. When it makes up a small portion of the total 
emission, the parameters B and Yradio can only be loosely 
constrained. The fitting results in the case where no pa- 
rameter of Model III is fixed are listed in Table A1. The 


maximum-likelihood estimates of Yradio are 0.92+0-32 and 


133593) for all RLQs and FSRQs, which are close to unity. 
The quoted error bars indicate lo uncertainties. The fitting 
results of Zamorani (1984) and Worrall et al. (1987) also 
suggest that there is a linear correlation between the radio 
and X-ray luminosities for the jet component. Miller et al. 
(2011) prefer a scenario where the excess X-ray emission 
of RLQs relative to RQQs is jet-linked, and this jet-linked 
X-ray emission is beamed with a smaller bulk Lorentz factor 
than that of the radio emission, which suggests Yradio < 1. 
However, FSRQs are more X-ray luminous than SSRQs at 
given Lyy and R (see Fig. 6), which is inconsistent with the 
idea that the jet-linked X-ray emission has less anisotropy 
than the radio emission. Note that fixing Yradio = 1 does not 
affect the result that the jet component is a very minor term 
but allows for an estimate of B with smaller uncertainty. 
The fitting for SSRQs in Table A1 is unphysical. Firstly, 
A € B would suggest that the jets are equal to or more im- 
portant than the corona in explaining the X-ray luminosities 
of SSRGs, which is in contrast with other samples. In fact, we 
expect much less contribution from the jets in SSRQs than 
in FSRQs. Secondly, the value of Tuy = yuv ^ Ya, = 1.01+9-15 
is not even close to y. Model II has similar issues when it 
is applied to SSRQs. The fitting results using Model II with 
all four free parameters are listed in Table A2. They are 
non-physical for the same reasons as the case of Model III. 
These issues with Model II and Model III (without fixing 
parameters) already suggest that they are not appropriate 
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models for SSRQs. To make a meaningful comparison be- 
tween a model that invokes a distinct jet component and 
the model attributing X-ray luminosity to the corona, Yuy of 
Model II is fixed to 0.63 in Table 4. 


This paper has been typeset from a TEX/IATEX file prepared by 
the author. 
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Table A1. The fitting results using Model III with all parameters allowed to vary. 


a = Y radio p Yuy Yradio 
Sample III: log Laxey = log (are Lo SQ0À + Bie) 
A B Yuv Yradio JY adio 
0.01 0.19 0.15 0.32 0.04 
All RLQs 0.75752) 0.03+9-19 05075 0.9293 01770 04 
FSRQs 0.777914 0.03756) 0.517903 1.267944 0.187) Qi 
SSRQs 0.32:0-3  0,35+0- 0.89*0.18 — 9,39*0-07 — 9 12:00 


—0.07 -0.11 —0.24 —0.06 


—0.08 


Table A2. Model fitting results for SSRQs using Model II with 
all parameters set free to vary. 


. = Yuv Yradio 
Sample II: log Loxey = log [An + Biles) 
A B Yuv Yradio 


0.07 0.03 0.09 0.05 
SSRQs 0.237544  0.39*9 5 — 1.037015 0.367 5753 
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